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CHAPTER  I 


INTRODUCTION 

Much  v$eful  information  concerning  the  properties  of  a  stellar 
atmosphere  can  be  obtained  from  a  study  of  the  intensities  of  the 
absorption  lines  appearing  in  the  stellar  spectrum,  A  commonly  employed 
method  of  analyzing  absorption  line  intensities  involves  the  application 
of  idealized  curves  of  growth.  The  purpose  of  this  paper  is  to  describe 
the  use  of  this  approach,  called  "grobanalyse"  by  German  workers,  in 
a  study  of  the  spectrum  of  the  star  0  Ursae  Majoris*. 

The  advantage  of  the  curve  of  growth  lies  in  the  fact  that  it 
may  be  used  to  handle  many  observational  data,  namely,  the  equivalent 
widths  of  the  lines.  The  equivalent  width  of  a  line  is  a  readily 
obtainable  quantity,  whereas  the  actual  shape  or  profile  pf  a  line, 
which  is  required  for  the  use  of  the  more  refined  model  atmosphere 
techniques,  is  difficult  to  determine,  especially  for  weak  lines.  The 
disadvantage  of  the  curve  of  growth  arises  from  its  being  a  statistical 
procedure  which  ignores  information  that  might  be  obtainable  from  the 
shapes  of  certain  lines. 

The  curve  of  growth  was  introduced  by  Minnaert  and  Slob  (1931)  and 
developed  and  applied  by  many  other  workers.  Menzel  (1936)  first 


*9  Ursae  Majoris  is  listed  by  Johnson  and  Morgan  (1953)  as  a 
subgiant  (spectral  type  F6,  luminosity  class  IV),  Its  visual  magnitude 
is  3.3  and  its  coordinates  are  Qr(1900)  m  9^26"',  6(1900)  ■  +52°8' 

(Keenan  and  Morgan,  1951)  • 
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developed  the  theory  as  it  applies  to  the  Schuster-Schwarzschild  (S-S) 
model  atmosphere'-f’wherein  the  entire  continuous  spectrum  is  produced 
by  a  sharply  defined  radiating  surface  or  photosphere,  overlaid  by  a 
reversing  layer  which  produces  all  the  line  absorption  but  no  continuous 
absorption.  Such  an  approximation  can  be  considered  quite  representative 
for  the  lines  of  neutral  metals  (e,g.,  Ca  I  and  Na  I  in  the  sun)  where 
the  number  of  absorbers  diminishes  with  optical  depth  because  of 
increased  ionization  with  increased  temperature.  The  other  extreme 
approximation,  where  the  lines  originate  in  the  same  layers  as  does  the 
continuous  spectrum,  is  embodied  in  the  Milne-Eddington  (M-E)  approxima¬ 
tion  of  a  ratio  of  line  fo  continuous  absorption  coefficient  which  is 
constant  with  depth.  Stromgren  (1937)  has  found,  at  least  in  some 
cases,  that  the  variation  of  this  ratio  with  optical  depth  is  sufficiently 
small  to  permit  a  mean  value  to  be  employed  over  the  layers  important 
in  line  formation,  in  which  case  the  M-'E  model  can  be  used.  The  true 
situation  for  any  line  lies  between  the  two  extremes  represented  by 
the  S-S  and  M-E  models. 

Wrubel  (1949,  1950,  1956)  has  computed  theoretical  curves  of  growth 
based  upon  the  M-E  model  for  pure  scattering  and  for  pure  absorption. 

The  processes  of  pure  scattering  and  pure  absorption  are  the  two  extreme 
mechanisms  by  which  a  line  may  be  formed.  Also,  Wrubel  (1954a)  has 
published  curves  for  the  S-S  model  with  pure  scattering.  All  of  Wrubel 's 
curves  of  growth  are  based  upon  Chandrasekhar's  exact  solutions  of  the 
equation  of  transfer  and  therefore  are  regarded  as  the  most  accurate 
available  for  the  above  mentioned  models.  Theoretical  curves  of  growth 
for  the  S-S  and  M-E  models  with  pure  absorption  have  been  calculated  by 
Hunger  (1956). 
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Although  these  Idealized  curves  can  be  brought  into  approximate 
Qoincidence  (Wrubel,  1954b,  and  Hunger,  1955),  significant  differences 
exist.  However,  the  scatter  in  the  observations  is  usually  large 
enough  that  it  is  difficult  to  determine  with  any  degree  of  certainty 
which  of  the  idealized  curves  is  best  fit  by  the  observations;  hence, 
it  was  deemed  edvisable  to  apply  all  four  sets  of  curves  to  the  study 
of  the  spectrum  of  9  Ursae  Majoris.  The  fqur  sets  of  curves  applied 
in  this  study  aye  thpse  for  pure  scattering  by  Wrubel  for  the  M-E  and 
S-'S  models  and  Hunger's  M’E  and  S- S  pure  absorption  curves.  For  reasons 
to  be  discussed  in  Chapter  II,  Hunger's  li»E  pure  absorption  curves 
wete  used  rathet  than  those  by  Wrubel  for  this  same  model. 

This  investigation  was  restricted  to  those  lines  for  which  labora¬ 
tory  absolute  fryalues  are  reasonably  well  known.  Recently,  Corliss  and 
Bozman  (1962)  published  a  list  of  experimentally  determined  absolute 
frvalues  for  25,000  spectral  lines  of  seventy  elements.  Due  to  the 
availability  of  this  extensive  tabulation,  there  appeared  to  be  no 
necessity  for  using  theoretically  calculated  f-values;  such  quantities 
determined  by  laboratory  measures  are  considered  to  be  more  reliable 
that  those  calculated  from  theory.  Other  experimentally  determined  f- 
values  employed  in  this  investigation  are  those  of  King  and  King  (1938) 
for  Fe  I  and  Ti  I,  Hill  and  King  (1951)  for  Cr  I,  and  Carter  (1949)  for 
Fe  1.  Also,  the  solar  log  of  Wright  (1948)  were  used  to  augment  the 

Ti  II  data. 

Only  lin^s  with  wavelengths  in  excess  of  4000  angstroms  were  used, 
since  the  spectrum  is  quite  complex  to  the  violet  of  A.4000.  Also,  many 
lines  in  this  region  are  on  the  broad  wings  of  the  H  and  K  lines  of 


ionized  calcitun, 
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Equivalent  width?  have  been  measured  for  359  lines  of  Ca  i,  Ti  I, 

Ti  II,  Cr  1,  Mn  I,  Fe.  I,  Co  I,  and  Ni  I.  The  spectral  region  covered 
is  WAOQO- 6500,  and  only  those  lines  which  were  sufficiently  well 
resolved  apd  unblended  were  used.  Kinetic  velocities,  damping  constants, 
excitation  temperatures,  and  abundances  were  determined  using  the 
thepretfcal  curves  of  growth  based  upon  both  the  M-E  and  S^S  models. 
Results  were  obtained  for  the  extreme  cases  of  pure  scattering  and 
pure  absorption  for  the  two  models. 


CHAPTER  XI 


THE  CURVE  OF  GROWTH 

Pef inition 

The  hurv^  of  growth  is  the  relation  between  the  intensity  of  an 
absorption  line  and  the  number  of  absorbing  atoms  active  in  producing 
the  Hne, 

The  Intensify  of  an  Absorption  Line 

The  term  ''intensity"  here  means  the  total  amount  of  energy  sub¬ 
tracted  from  the  continuous  spectrum  by  the  absorption  line  and  is  best 
specified  by  a  quantity  called  the  "equivalent  width,"  W,  of  the  line. 
The  profile  of  a  line  is  the  plot  of  the  flux  at  each  point  in  the 
line  versus  the  wavelength.  Hence,  the  area  enclosed  by  the  profile 
and  the  continuum  is  the  total  amount  of  energy  absorbed  in  the  line. 
The  equivalent  width  is  the  width  of  the  rectangle  having  this  same 
area  and  with  a  height  equal  to  that  of  the  adjacent  continuum.  These 
relations  are  shown  in  Figure  1. 

The  advantages  of  expressing  the  intensity  of  a  line  in  terms  of 
the  equivalent  width  are  that  W  is  much  less  affected  by  the  finite 
resolution  of  the  spectrograph  than  is  the  line  profile,  and  W  is 
in4ependent  of  the  Doppler  effects  of  stellar  rotation  and  large  scale 
turbulence  in  the  stellar  atmosphere. 
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Figure  1.  The  Equivalent  Width,  W,  of  a  Spectral  Line. 


The  Number  of  Absorbing  Atoms 


The  number  of  absorbing  atoms  for  a  line  originating  from  a  level 
n  and  tenninating  on  level  j  is  given  by  where  is  the  number 

of  atoms  in  level  n,  and  is  the  f-value  or  oscillator  strength  of 
the  transition  from  level  n  to  level  j.  When  an  atom  in  level  n  is 
ejcposed  to  radiation  of  all  frequencies  it  may  absorb  any  one  of 
several  diiierent  quanta.  In  specifying  the  number  of  absorbing  atoms 
for  a  particular  line  it  is  necessary  to  include  the  f-value  to  take 
into  account  the  probability  of  the  atom  absorbing  the  particular 
quantum  corresponding  to  the  transition  in  question.  The  f- values 


depend  only  upon  the  structure  of  the  atom,  and  may  be  computed  from 
theory  or  measured  in  the  laboratory. 


Observed  Curves  of  Growth 

Lines  Produced  in  an  Absorption  Tube 

Let  us  consider  the  growth  of  a  particular  line  as  the  number  of 
absorbers  is  increased.  This  "growth"  of  a  line  may  be  observed  in  the 
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laboratory  by  varying  tha  concentration  of  the  element  under  investiga¬ 
tion  in  an  absorption  tube.  Lines  with  strengths  depending  upon  the 
concentration  ate  formed  when  light  from  an  incandescent  source  shines 
through  the  tube.  When  the  number  of  absorbers  is  small  the  principal 
contributor  to  the  equivalent  width  is  thermal  Doppler  broadening-,  with 
the  result  that  W  is  proportional  to  Nf.  Hence,  plotting  log  W  ("log" 
represents  the  logarithm  to  the  base  ten)  versus  log  Nf  for  small 
values  of  Nf  gives  a  portion  of  the  curve  of  growth  which  is  linear 
and  has  slope  equal  to  unity.  Increasing  the  number  of  absorbers  causes 
the  center  of  the  line  to  become  saturated,  and  a  linear  relation ; between 
W  and  Nf  no  longer  exists.  The  point  of  departure  from  linearity  is 
called  the  "Knee"  of  the  curve  of  growth.  Beyond  the  knee  W  increases 
slowly,  being  proportional  to  the  square  root  of  the  natural  logarithm 
of  Nf .  As  the  number  of  absorbers  is  increased  still  more,  the-  damping 

portion  of  the  absorption  coefficient  is  reached.  Here  the  number  of 

atoms  is  large  enough  to  sufficiently  perturb  the  energy  levels  so  that 

an  appreciable  amount  of  absorption  occurs  away  from  the  saturated 

center  of  the  line.  Thus -the  line  begins  to  grow  more  rapidly.  In 
this  portion  of  the  curve  of  growth  W  is  proportional  to  the  square 
root  of  Nf. 

This  discussion  assumes  that  the  temperature  of  the  absorbing  atoms 
is  a  particular  constant  value.  The  position  of  the  knee  of  the  curve 
of  growth  varies  with  the  temperature.  As  the  temperature  increases, 
the  Maxwellian  velocity  distribution  becomes  more  extended,  and  satura¬ 
tion.  occurs  at  a  greater  line  width,  giving  a  higher  knee  on  the  curve 
of  growth.  Also,  the  point  at  which  the  damping  portion  of  the  curve 
begins  moves  closer  to  the  knee  as  the  vi  l-’e  of  the  damping  parameter 
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increases.  The  demping  parameter,  a,  is  the  ratio  of  the  effective 
natural  line  width  to  the  Doppler  width.  Figure  2  illustrates  these 
relations , 


(a)  (b) 

Figure  2.  Effects  of  Temperature  (a)  and  Damping  (b)  on  the  Curve  of 
Growth. 


Lines  Produced  in  a  Stellar  Atmosphere 

The  discussion  thus  far  has  dealt  with  the  growth  of  a  single 

line  as  the  number  of  absorbing  atoms  is  increased.  A  single  line 

cannot  be  observed  to  grow  in  a  stellar  spectrum  because  does 

not  change.  However,  in  going  from  one  line  to  another  of  the  same 

element.  N  f  .  does  change,  since  f  .  varies  from  line  to  line,  and  N 
’  n  nj  - -  ®  ’  nj  ’  n 

also  varies  if  the  lines  arise  from  different  levels.  Again  a  curve 

of  growth  results  from  the  relation  between  the  equivalent  widths  of 

the  various  lines  and  the  corresponding  values  of  N  f  .. 

n  nj 

In  beginning  a  study  of  a  stellar  spectrum,  the  values  of  are 
not  known,  and  one  of  the  objectives  of  a  curve  of  growth  study  is  to 
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determine  these  quantities.  The  procedures  employed  for  this  purpose 
are  described  later  in  this  chapter. 

In  applying  the  curve  of  growth  technique  to  many  different  lines 
of  an  element  it  is  necessary  to  assume  that  all  the  lines  are  formed 
at  the  same  temperature.  The  temperature  in  a  stellar  atmosphere 
increases  roughly  as  T  ,  where  T  is  the  optical  depth.  Thus  the 
assumption  is  made  that  all  the  lines  are  formed  at  a  particular  depth 
in  the  atmosphere.  This  may  be  a  poor  assumption  for  weak  lines  and 
the  wings  of  strong  lines  where  the  absorption  coefficient  is  small. 

In  this  case  the  line  is  formed  in  a  thick  layer  through  which  the 
temperature  may  change  appreciably.  In  the  center  of  a  strong  line 
the  absorption  coefficient  is  large,  and  the  line  is  effectively  formed 
in  a  thin  l^yer  in  which  the  temperature  may  be  practically  constant. 
Considering  non- resonant  lines  of  a  neutral  atom,  the  temperature  in 
the  outer  layers  of  the  atmosphere  may  be  too  low  to  appreciably  excite 
the  atoms  to  levels  capable  of  absorbing  these  lines.  In  very  deep 
layers  of  the  atmosphere  ionization  will  reduce  the  number  of  absorbers 
to  a  negligible  number.  The  lines  are  consequently  formed  in  inter¬ 
mediate  layers,  and  the  temperature  in  this  case  is  characteristic  of 
these  layers. 

It  is  also  necessary  to  assume  that  a  particular  value  of  the 
damping  parameter  can  be  taken  as  representative  for  all  the  strong 
lines.  This  assumption  is  somewhat  questionable,  since  the  damping 
parameter  is  known  to  differ  from  line  to  line.  However,  this  effect 
does  not  come  into  play  for  weak  lines,  and  it  is  usually  possible  to 
determine  a  value  for  the  damping  parameter  that  fits  the  observations 


reasonably  well . 
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As  mentioned  earlier,  the  position  of  the  knee  of  the  curve  of 
growth  depends  upon  the  temperature.  Another  factor  which  may  influence 
the  position  of  the  knee  is  microturbulence,  i.e.,  turbulent  motions  of 
the  gases  of  the  atmosphere  where  these  motions  or  eddies  are  confined 
within  the  layer  in  which  the  lines  are  produced.  The  presence  of  micro¬ 
turbulence  increases  the  widths  of  the  lines  thus  increasing  the  number 
of  absorbers  required  for  saturation  of  the  line  center.  Hence,  the 
position  of  the  knee  is  elevated.  Effects  of  microturbulence  are 
especially  evident  in  curves  of  growth  for  giant  and  supergiant  stars. 
Struve  and  Elvey  (1934)  first  called  attention  to  this  phenomenon. 

Theoretical  Curves  of  Growth 

i'  " 

Calculations  of  theoretical  curves  of  growth  are  based  upon  various 
assumptions  as  to  the  structure  of  stellar  atmospheres  and  the  mechanism 
of  line  formation. 

The  Schuster- Schwarzschild  Model 

The  Schuster- Schwarzschild  (S-S)  model  depicts  the  atmosphere  as 
consisting  of  a  "reversing  layer"  of  depth  H  overlying  a  sharply  defined 
radiating  surface  or  photosphere.  The  continuous  spectrum  is  produced 
below  the  photosphere,  and  the  reversing  layer  is  the  region  of  line 
formation.  The  optical  depth  of  the  reversing  layer  is  a  function  of 
the  absorption  coefficient  of  the  line,  and  there  is  no  continuous 
absorption  in  the  reversing  layer.  A  portion  of  the  radiation  incident 
upon  the  photosphere  from  below  is  transmitted  while  the  remainder  is 
diffusely  reflected.  The  transmitted  energy  is  subject  to  selective 
absorption  and  scattering.  Radiation  returned  to  the  photosphere  is 
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absorbed  there  by  continuous  absorption.  For  a  point  in  the  wing  of  a 
line  the  optical  thickness  of  the  reversing  layer  is  small,  while  it  is 
much  greater  in  the  line  core.  This  model  is  a  good  approximation  for 
neutral  metals  where  the  number  of  absorbers  diminishes  rapidly  with 
depth  because  of  increased  ionization  with  increased  temperature. 

The  Milne-Eddington  Model 

In  the  Milne-Eddington  (M-E)  model  the  lines  and  the  continuum  are 

formed  in  the  same  region.  For  each  frequency  the  ratio  of  the  line 

absorption  coefficient  to  the  continuous  absorption  coefficient  (J?  /k  = 

V  V 

71  )  is  assumed  to  be  constant  with  depth,  although  H  and  k  may  indi- 

'y  V  V 

vidually  change  with  depth.  In  the  line  wings  there  are  few  absorptions 
and  reemissions,  and  consequently  the  emergent  radiation  follows  a 
short  path  through  the  atmosphere.  In  the  line  center  the  quanta 
follow  longer  paths  (many  absorptions  and  reemissions)  and  the  probabil¬ 
ity  of  continuous  absorption  is  greater  (say  due  to  H  in  stars  like 
the  sun)  and  the  energy  may  be  lost. 

Thus,  in  both  models  a  selective  process  followed  by  continuous 
absorption  forms  the  line.  Figure  3  illustrates  these  two  idealized 
models . 

: . 

. .4 


S-S  Model  M-E  Model 

Figure  3.  The  S-S  and  M-E  Models.  Dots  represent  selective  absorption 
and  emission;  hatches  represent  continuous  absorption. 
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The  Formation  of  a  Line 

The  mechanism  of  spectral  line  formation  is  determined  by  conditions 
existing  in  the  atmosphere.  If  the  atmosphere  is  in  a  state  of  radiative 
equilibrium,  then  for  each  quantum  absorbed  a  like  quantum  is  emitted. 

This  process  is  called  "scattering."  Scattering  is  coherent  if  the 
frequency  of  the  emitted  quantum  is  exactly  the  same  as  the  frequency  of 
the  absorbed  quantum.  If  the  frequency  of  the  emitted  quantum  is  slightly 
different  from  that  of  the  absorbed  quantum,  the  scattering  is  non¬ 
coherent.  If  the  atmosphere  is  in  "local  thermodynamic  equilibrium," 
i.e.,  Kirchhoff's  law  holds  at  each  point  of  the  atmosphere  for  the 
local  temperature,  then  the  lines  are  formed  by  the  process  of  "absorp¬ 
tion."  Here  there  is  no  one-to-one  correspondence  between  the  absorbed 
frequency  and  the  emitted  frequency,  and  the  absorbed  quantum  is  "lost." 

Hence,  a  spectral  line  may  be  formed  by  scattering  or  absorption 
or  by  some  combination  of  the  two, 

The  Theoretical  Curves  Used  in  This  Study 

The  four  sets  of  theoretical  curves  used  in  this  investigation  are 
as  follows:  1)  M-E  pure  scattering  by  Wrubel,  2)  S-S  pure  scattering  by 
Wrubel,  3)  M-E  pure  absorption  by  Hunger,  and  4)  S-S  pure  absorption  by 
Hunger . 

Curves  for  M-E  pure  absorption  have  also  been  calculated  by  Wrubel 
and  published  in  an  article  by  Aller  (1960) .  The  calculations  were 
based  upon  Chandrasekhar’s  exact  solution  of  the  equation  of  transfer 
and,  therefore,  should  be  the  most  accurate  available  for  the  M-E  pure 
absorption  case.  However,  these  curves  include  only  two  values  (-1  and 
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-3)  of  log  a,  where  a  is  the  damping  parameter.  Because  of  this  limita¬ 
tion  they  have  not  been  employed  in  this  study. 

The  integral  equation  for  the  equivalent  width,  W,  is 

+00 

^  =  Jd  -  R)du-,  (1) 


where  b  is  the  Doppler  width  of  the  line,  R,  the  "residual  intensity,” 
is  the  ratio  of  the  flux  of  the  radiation  at  a  point  in  the  line  to  the 
flux  in  the  adjacent  continuum,  and  w  is  the  distance  from  the  center  of 
the  line  in  units  of  the  Doppler  width.  Here  u"  =  (X  -  X^)/b,  where  X 
is  the  wavelength  at  a  point  in  the  line,  and  X^  is  the  wavelength  at 
the  center  of  the  line.  Also,  b  =  A.^v/c,  where  v  is  the  most  probable 
velocity  of  the  atoms  and  c  is  the  velocity  of  light. 

It  should  be  noted  that  the  velocity  v  includes  both  the  effects 
of  thermal  and  microturbulent  motions.  That  is. 


'i^h  \urb« 


where  v^^  is  the  most  probable  thermal  velocity,  and 


(2) 

is  the  most 


probable  velocity  of  the  microturbulence.  It  is  assumed  that  the  turbur 
lent  elements  have  a  random  distribution  of  velocities.  The  thermal 
velocity  is  given  by 

\h  'Nl  M  ’  (3) 


where  k  is  Boltzmann's  constant,  T  is  the  kinetic  temperature,  and  M 
is  the  mass  of  the  atom  producing  the  lines. 

The  residual  intensity,  R,  is  found  from  the  solution  of  the  equation 
describing  the  transfer  of  radiation  through  a  stellar  atmosphere.  R 


depends  upon  assumptions  as  to  the  atmospheric  model  and  the  mechanism 
of  line  formation.  In  his  calculations  of  theoretical  curves  of  growth. 
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Wrubel  employed  the  expressions  for  R  resulting  from  Chandrasekhar's 
exact  solutions  of  the  equation  of  transfer.  Expressions  for  R  which 
follow  from  the  work  of  Unsold  (1955)  were  used  by  Hunger  in  his  calcu¬ 
lations  of  theoretical  curves  of  growth. 

Both  Wrubel 's  and  Hunger's  curves  employ  the  assumption  that  the 
Planck  function  varies  linearly  with  optical  depth  in  the  continuum; 
i.e.,  By  (T)  =  B°  +  B^Ty.  This  leads  to  a  limb  darkening  law  in  the 
continuum  that  is  linear  in  cos  9,  where  0  is  the  angle  between  the  line 
of  sight  and  the  normal  to  the  stellar  surface.  In  the  S-S  model  this 
is  expressed  as  I  =  1°  +  cos  9,  and  in  the  M-E  model  as  I  =  B°  +  B^cos  0. 
The  effect  of  the  addition  of  an  exponential  term  to  the  linear  expression 
for  the  Planck  function  has  been  investigated  by  Hunger  (1956)  .  He  con¬ 
cluded  that  the  addition  of  such  a  term  has  negligible  effect  and  is 
therefore  unnecessary. 

The  expression  for  the  residual  intensity,  as  given  by  Chandrasekhar 
and  used  in  Wrubel 's  calculations,  involves  the  ratio  of  the  limb  darken¬ 
ing  coefficients,  I°/I^  or  B°/B^.  Pure  scattering  curves  are  available 
for  the  S-S  model  with  I°/I^  equal  to  h,  1,  and  2,  and  for  the  M-E 
model  with  B°/B^  equal  to  1/3,  2/3,  4/3,  and  10/3. 

Equation  (1),  with  the  appropriate  expressions  for  R,  was  employed 
by  Wrubel  in  his  calculations.  Consequently,  the  ordinate  of  his  curves 
is  log  W/b.  Hunger  divided  Equation  (1)  by  R^,  the  "limiting  depth"  of 
the  line,  to  obtain 


where 


_ W_ 

2R  b 
c 


(4) 
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R  = 
c 


3  ’ 


(5) 


1  +  t  ^ 
2 


and  the  symmetry  of  the  integral  about  kt  equal  to  zero  has  been  taken 

into  account.  The  ordinate  of  Hunger's  curves  is  log  W/2R  b.  For  the 

c 

cases  of  M-E  and  S-S  pure  absorption,  the  ratio  (1-R)/R^  is  independent 

of  the  limb  darkening  coefficients.  This  results  from  the  particular 

expressions  for  R  used  by  Hunger.  Therefore,  Hunger's  formulation  of 

the  curve  of  growth  for  M-E  and  S-S  pure  absorption  avoids  the  necessity 

o  1 

of  having  to  calculate  separate  curves  for  each  value  of  B  /B  .  It  is 
to  be  noted,  however,  that  the  limb  darkening  is  still  taken  into  account 
by  Hunger  since  B°/B^  is  involved  in  the  ordinate  of  his  curves. 

In  the  calculations  of  theoretical  curves  of  growth  the  damping 
constant,  a,  is  considered  to  be  a  free  parameter.  Wrubel's  pure  scat¬ 
tering  curves  were  calculated  for  values  of  log  a  running  from  -3.0  to 
-1.0  in  steps  of  0.4.  The  pure  absorption  curves  published  by  Hunger 
were  computed  for  values  of  log  a  running  from  -4.3  to  +0.7  in  steps  of 
0.5.  He  also  included  a  curve  for  log  a  =  -<». 

For  the  M-E  model  (both  scattering  and  absorption)  R  is  a  function 
of  Tlp>  where 


L 


^O  H 


(6) 


H  is  the  fictitious  absorption  coefficient  at  the  center  of  the  line, 
o 

and  K  is  the  mean  continuous  absorption  coefficient.  The  abscissa  for 
both  the  M-E  scattering  and  M-E  absorption  curves  is  log  T]^.  In  Hunger's 
notation  the  abscissa  is  log  C,  but,  as  he  point-s.  out,  C  =  for  M-E 
pure  absorption. 

For  the  transition  from  level  n  to  level  j,  is  given  by 
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u  TT  e  A  N  f  , 

-n  =  _ _ n_ni 

°  H  2,  _ 

me  bp  K 


where 


e  =  charge  of  the  electron, 

m  =  mass  of  the  electron, 

A  =  wavelength  of  the  line  in  question, 

N  =  number  of  atoms  in  level  n  per  unit  volume, 
n 

f  .  =  oscillator  strength  of  the  transition  from 
nj 

level  n  to  level  j, 
c  =  speed  of  light, 
b  =  Doppler  width  of  the  line,  and 
p  =  density  of  the  stellar  material. 

Here  Tt  is  the  mean  absorption  coefficient  per  gram  of  stellar  material. 
Using  the  relation  b  =  Av/c,  it  follows  that 


log  Tl  =  log  N  / pH  +  log  C  +  log  f  .A  +  log  c/v, 
'o  nj  0  7 


where 


The  residual  intensity,  R,  in  the  S-S  model  is  a  function  of  t  the 

o 

fictitious  optical  depth  at  the  center  of  the  line.  For  S-S  pure  scatter¬ 


ing  Wrubel  uses  log  T  as  abscissa  while  Hunger,  for  S-S  pure  absorption, 

3 

uses  log  C,  where. C  =  in  this  case.  Hence,  Hunger's  abscissa  differs 

from  that  of  Wrubel  by  the  additive  constant  log  3/2. 

^22 

_  ^  %  ^XT  £  Tl 


T  =  £  pH  = 
o  *o^ 


TT^e  A  N  f  ,H 
n  n  1 


where  H  is  the  depth  of  the  reversing  layer.  Then 


log  *  log  N^H  +  log  C  +  log  c/v  . 


(11) 
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The  Idealized  Situation 

If  1)  the  atmosphere  were  built  on  either  the  M-E  or  the  S-S  model, 

2)  the  temperature  and  density  were  constant  throughout  the 
atmosphere , 

3)  turbulence  were  negligible, 

4)  the  effects  of  broadening  were  constant  with  depth, 

5)  measurements  for  many  lines  from  a  single  level  were  available, 

6)  the  lines  were  formed  either  by  pure  scattering  or  pure 
absorption,  end 

7)  absolute  f-values  were  known, 

then  1)  the  damping  constant  could  be  found, 

2)  the  temperature  could  be  determined,  and 

3)  the  abundance  of  the  element  could  be  determined, 
all  from  the  curve  of  growth. 

For  such  an  idealized  situation  the  procedure,  for  the  S-S  pure 

scattering  model  for  example,  would  be  as  follows:  plot  log  W/X  against 

log  fX  and  superpose  this  observed  curve  of  growth  upon  the  theoretical 

curves  drawn  to  the  same  scale;  shift  vertically  and  horizontally  with 

no  rotation  until  the  curves  are  brought  into  coincidence;  obtain  log  a 

by  noting  which  theoretical  curve  best  fits  the  damping  portion  of  the 

observed  curve;  from  the  vertical  shift  determine  log  c/v  and  thus  v; 

from  this  value  of  v  calculate  the  kinetic  temperature  using  Equation  (3), 

since  v  =  v^,  in  the  absence  of  turbulence;  from  the  horizontal  shift, 
th 

which  gives  log  N^H  +  log  C  +  log  c/v,  find  log 

Departures  from  the  idealized  conditions  may  place  severe  limitations 
on  the  curve  of  growth  technique. 
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Reduction  of  Curve  of  Growth  Data 


The  Boltzmann  Equation 


In  practice,  a  large  number  of  lines  all  arising  from  the  same 
level  of  excitation  usually  are  not  available.  Hence,  lines  arising 
from  different  levels  are  plotted  on  the  same  curve  of  growth.  For 
this  purpose  the  assumption  is  made  that  an  excitation  temperature  can 
be  specified  so  that  the  distribution  of  atoms  in  the  various  energy 
levels  is  given  by  the  Boltzmann  equation,  i.e., 

(12) 


Ngn 

\=  —  exp(-Xg/kT), 


where 


N  =  total  number  of  atoms  per  unit  volume  of  the  element 
in  question, 

g^  =  2J  +  1  =  statistical  weight  of  level  n, 

V  =  excitation  potential  of  level  n, 
k  =  Boltzmann's  constant, 

T  =  excitation  temperature,  and 

u  =  the  partition  function  =  )g.w  exp(-x  . /kT)  summed 
over  all  levels  in  the  atom,  w  is  a  correction  for 
pressure  effects. 

By  replacing  in  Equations  (8)  and  (11)  by  the  expression  in 
Equation  (12),  it  follows  that 


and 


log  71q  =  log  N/pn  +  log  C  +  log  "  log  u  +  log  c/v  -  0Xg>  (13) 


log  =  log  NH  +  log  C  +  log  "  log  o  +  log  c/v  -  9Xg. 


where  9  ■  5040/T,  if  Xg  Is  expressed  in  electron  volts, 


(14) 
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The  Method  of  Analysis 

For  purposes  of  the  following  discussion,  let  us  confine  our 
attention  to  the  application  of  the  S-S  pure  scattering  model  to  the 
analysis  of  the  lines  belonging  to  a  single  element  in  a  particular 
stage  of  ionization. 

Notice  that  all  the  terms  on  the  right-hand  side  of  Equation  (14), 
except  log  gfX*  and  constant  for  all  the  lines.  For  a  group  of 

lines  that  all  arise  from  the  same  atomic  term,  is  practically  con¬ 
stant;  therefore,  log  gfX  is  the  only  variable  quantity  involved  in  the 
abscissa  of  a  curve  of  growth  limited  to  such  a  group.  The  ordinate  of 
the  theoretical  curves  can  be  written  as 

log  W/b  =  log  W/X  +  log  c/v,  (15) 

where  log  c/v  is  assumed  to  be  constant  for  all  the  lines.  Hence,  plots 
that  can  be  compared  to  the  theoretical  curves  can  be  obtained  as  follows 
all  lines  arising  from  the  same  atomic  term  are  plotted  together  on 
transparent  paper  with  log  gfX  as  abscissa  and  log  W/X  as  ordinate. 

Each  plot  obtained  in  this  manner,  i.e.,  by  taking  the  lines  in 
groups  corresponding  to  particular  atomic  terms,  defines  a  portion  of  the 
observed  curve  of  growth.  These  plots  are  then  superposed  upon  plots  of 
the  theoretical  curves,  drawn  to  the  same  scale,  and  shifted  vertically 
and  horizontally  (with  no  axial  rotation)  until  the  best  fit  is  obtained. 
If  enough,  strong  lines  are  available,  a  value  of  log  a  can  be  determined 
by  choosing,  from  among  the  theoretical  curves  for  various  values  of 
log  a,  the  curve  that  best  fits  the  plots  containing  the  strong  lines. 


*To  simplify  the  notation,  the  subscripts  on  g  and  f  will  be  omitted 
henceforth . 
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If  plots  are  available  which  fit  around  the  knee  of  the  curve  of  growth, 
it  is  possible  to  obtain  log  c/v  from  the  vertical  shift--since 

vertical  shift  =  log  W/b  -  log  W/X  =  log  c/v.  (16) 

It  should  be  noted  that  the  same  vertical  shift  must  be  used  for  all 
plots,  since  log  c/v  is  assumed  to  be  the  same  for  all  the  lines. 

The  horizontal  shift  yields  Alog  X,  where 

Alog  X  =  log  -  log  gfX.  (17) 

Also, 

Alog  X  =  [log  NH  +  log  C  -  log  u  +  log  c/v]  -  9Xg> 
which  can  be  seen  by  referring  to  Equation  (14).  In  Equation  (18),  the 
quantity  in  brackets  is  constant  for  all  the  lines.  Hence,  plotting 
Alog  X  versus  Xg  each  of  the  different  terms  should  give  a  straight 
line  with  slope  equal  to  The  excitation  temperature  can  be  found 

since  T  =  5040/9 . 

From  Equation  (18)  it  is  apparent  that  Alog  X  +  Qx^  is  equal  to  a 
constant,  henceforth  called  the  "shift."  Then 

shift  =  log  NH  +  log  C  log  u  +  log  c/v.  (19) 

Consequently,  the  abundance  is  given  by 

log  NH  =  shift  -  log  C  +  log  u  -  log  c/v.  (20) 

From  the  values  of  the  shift  and  log  c/v  determined  as  described  above, 
log  NH  can  be  calculated  from  Equation  (20)  with  the  aid  of  log  C, 
which  is  a  known  constant  (-12.30),  and  log  u,  which  has  been  tabulated 
by  Aller  (1960) .  Since  N  is  the  number  of  atoms  per  cubic  centimeter 
and  H  is  the  depth  of  the  reversing  layer,  NH  is  the  number  of  atoms  in 
a  one  square  centimeter  column  of  the  reversing  layer. 

The  above  discussion  describing  the  application  of  the  S-S  pure 
scattering  model  can  be  made  valid  for  the  other  three  models  by  a  few 
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simple  alterations.  These  changes  are  outlined  in  the  paragraphs  to 
follow. 

The  discussion  relating  to  the  application  of  the  S-S  pure  scatter¬ 
ing  model  can  be  made  valid  for  M-E  pure  scattering  by  replacing  log  gfX 
by  log  gfX  +  Alog  by  T]^,  and  NH  by  N/pn.  The  addition  of 

Alog  T|^  to  log  gfX  takes  into  account  limb  darkening  effects.  As 
mentioned  earlier,  Wrubel's  curves  have  been  calculated  for  several 
values  of  or  with  a  separate  curve  resulting  for  each 

value.  Since  these  quantities  depend  upon  the  wavelength,  and  therefore 
vary  from  line  to  line,  each  line  should  fit  on  a  different  curve  of 
growcn.  This  effect  is  negligibly  small  for  the  S-S  pure  scattering 

model  and  need  not  be  considered.  For  M-E  pure  scattering,  however, 

0,1 

in  order  to  plot  all  lines  on  the  same  curve  of  growth,  say  for  B  /B  = 
2/3,  it  is  necessary  to  add  a  quantity,  Alog  Tj^,  to  the  value  of  log  gf\ 
for  each  line.  This  shifts  the  abscissa  of  the  line  by  the  proper  amount 
to  take  into  account  the  limb  darkening  effect.  These -limb  darkening 
corrections  can  be  determined  by  interpolating  in  the  tables  published 
by  Wrubel  for  the  case  of  M-E  pure  scattering. 

Hunger's  curves  for  M-E  and  S-S  pure  absorption  do  not  depend  upon 
the  limb  darkening  coefficients.  Limb  darkening  effects  are  taken  into 
account  by  the  inclusion  of  R  in  the  ordinates.  R  can  be  calculated 
for  each  line  from  Equation  (5),  using  appropriate  values  of  B°/B^. 

To  convert  the  discussion  relating  to  S-S  pure  scattering  to  one  apply¬ 
ing  to  M-E  and  S-S  pure  absorption  proceed  as  follows:  M-E  model-- 
replace  W  by  W/2R^,  by  T)^,  and  NH  by  N/pTt;  S-S  model--replace  W  by 

W/2R  ,  T  by  3t  /2,  and  "shift"  by  "shift"  +  log  3/2.  The  log  3/2 
CO  o 

comes  about  because  the  abscissa  here  is  log  3t  /2,  not  simply  log  T  . 

0  o 


22 


In  the  M-E  model  the  abundance  of  the  element  is  expressed  as 
N/pTi.  N/p  is  the  number  of  atoms  of  the  element  in  question  per  gram 
of  stellar  material,  and  1/h  is  a  measure  of  the  number  of  grams  of 
stellar  material  in  a  one  square  centimeter  column  of  the  atmosphere. 
Hence,  N/pK  is  a  measure  of  the  number  of  atoms  in  a  one  square  centi¬ 
meter  column  of  the  atmosphere. 

If  the  element  being  studied  exists  in  several  stages  of  ionization, 
then  the. total  abundance  is  given  by  the  Saha  equation.  In  logarithmic 
form, 

log  N^/N^  =  -  log  -  0x^  +  2.5  log  T  -  0,48  +  log  2uj^/u^,  (21) 

vhere  is  ionization  potential  in  electron  volts,  P^  is  the  elec- 

2 

tron  pressure  in  dynes  per  cm  ,  is  the  number  of  singly  ionized 

3  3 

atoms  per  cm  j  is  the  number  of  neutral  atoms  per  cm  ,  u^  is  the 

partition  function  of  the  singly  ionized  atoms,  and  u^  is  the  partition 

function  of  the  neutral  atoms.  The  equation  may  be  successively  applied 

to  include  all  significant  stages  of  ionization. 

Summary  of  Steps  Involved  in  a  Curve  of  Growth  Analysis 

1)  Obtain  spectrograms  of  the  star  to  be  investigated  .--It  is 
necessary  to  use  spectra  of  the  highest  dispersion  possible  so  that 
effects  of  blending  will  be  minimized. 

2)  From  these  plates  obtain  microphotometer  tracings .--Micro¬ 
photometer  tracings  are  plots  of  the  density  of  the  spectral  image 
versus  the  wavelength. 

3)  Obtain  intensitometer  tracings . '-The  intensitometer  converts 
the  density  of  the  image  into  intensity..  Calibration  involves  the  use 

of  sources  of  known  intensity.  The  central  portion  of  the  density-intensity 
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curve  for  the  plate  is  used. 

At  the  Dominion  Astrophysical  Observatory  steps  two  and  three  were 
performed  separately.  At  other  observatories,  Mt .  Wilson  for  example, 
the  two  may  be  combined  into  a  single  operation.  Line  profiles  may 
also  be  obtained  by  photoelectric  spectrophotometry.  The  instruments 
and  procedures  used  in  these  various  techniques  have  been  described  by 
Wright  (1962)  . 

4)  Draw  in  the  continuum. --The  location  of  the  continuum  is  a 
great  source  of  error.  The  continuum  can  be  drawn  directly  on  the 
intensitometer  tracing  or  first  drawn  on  the  microphotometer  tracing 
and  then  transferred.  The  procedure  employed  depends  upon  which  seems 
to  be  the  most  reliable  for  the  wavelength  region  concerned, 

5)  Identify  the  lines  and  approximate  the  prof iles .--If  the 
apparent  line  profile  is  irregular  or  asymmetric,  several  tracings  may 
be  compared  to  determine  if  the  variations  are  real,  A  possible  cause 
for  irregularities  and  asymmetry  is  blending.  In  drawing  the  profiles 
these  possible  effects  of  blending  should  be  taken  into  account.  Con¬ 
cerning  the  approximation  of  the  profiles  a  question  arises:  should 

the  profiles  be  drawn  so  that  the  wings  are  inversely  proportional  to  • 
2 

tf  ,  as  theory  indicates,  or  should  the  profiles  be  assumed  to  be 
triangular  in  shape  as  some  workers  do?  The  procedure  indicated  by 
theory  was  followed  in  this  investigation. 

6)  Determine  equivalent  widths. --The  following  steps  are  needed: 
a)  planimeter  the  lines  to  determine  the  area  enclosed  by  the  profile 
and  the  continuum,  b)  measure  the  height  of  the  continuum,  c)  measure 
the  dispersion,  and  d)  calculate  the  equivalent  width.  Deutsch  (1954) 
has  uiscussed  the  effect  of  the  finite  width  of  the  microphotometer 
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analyzing  slit  on  the  equivalent  width.  He  found  that  as  the  slit 
opening  increased  so  did  the  measured  equivalent  width. 

7)  Bring  together  lines  arising  from  the  same  term  and  plot 
log  W/\  (or  log  W/2R^X)  versus  log  gfX  on  transparent  paper .--Appro¬ 
priate  limb  darkening  corrections  should  be  included. 

8)  Fit  these  plots  to  the  theoretical  curves  of  growth, --Then, 
a)  the  theoretical  curve  which  is  best  fit  by  the  strong  lines  gives 
log  a,  b)  log  c/v  is  obtained  from  the  vertical  shift,  c)  the  excitation 
temperature  is  found  from  the  slope,  determined  by  the  least  squares 
method,  of  the  straight  line  that  best  fits  the  plot  of  Alog  X  versus 
y.g,  and  d)  the  abundance  is  determined  from  the  horizontal  shift. 

Steps  7  and  8  must  be  repeated  for  each  element  and  each  degree 
of  ionization. 

9)  Apply  the  Saha  equation  to  get  the  total  abundance  of  all 


observed  elements . 


CHAPTER  III 


OBSERVATIONAL  MATERIAL 

Spectrograms  and  Tracings 

The  9  Ursae  Majoris  spectrograms  used  in  this  study  were  taken  at 
the  Cassegrain  focus  of  the  72- inch  telescope  of  the  Dominion  Astro- 
physical  Observatory  by  Dr.  K,  0.  Wright,  Some  representative  spectro¬ 
grams  are  shown  in  Plate  I. 

The  dispersion  varied,  depending  upon  the  spectrograph  used.  Fcr 
the  Littrow  spectrograph  with  the  Wood  grating  (15,000  lines/inch) 
the  dispersion  was  approximately  7.5  A/mm  for  the  second  order  spectra 
in  the  range  \X4800-6750;  for  third  order  spectra  in  the  range  \X3750- 
4500,  the  dispersion  was  about  4,5  A/mm.  When  the  Bausch  and  Lomb 
grating  No.  496  (30,000  lines/inch)  was  used  in  the  second  order,  the 
dispersion  was  about  3.2  A/mm.  For  the  three-prism  spectrograph  the 
dispersion  varied  from  about  5  A/mm  to  15  A/rnm  over  the  wavelength  range 
studied . 

The  microphotometer  and  intensitometer  tracings  were  made  at 
Victoria  by  Dr.  Leon  W.  Schroeder.  The  magnification  of  the  tracings 
is  200.  Details  of  the  spectrograms  and  tracings  are  given  in  Table  I. 

Figure  4  is  a  reproduction  of  a  portion  of  one  of  the  microphotom¬ 
eter  tracings  showing  the  estimated  position  of  the  continuum  and  the 
profiles  of  some  representative  lines. 
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Plate 
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Spectrum  of  9  Ursae  Majoris  in  the  Regions  AX.4000-4100  and  4915-5260.  Wavelengths  and 
identifications  have  been  marked  for  representative  lines.  The  spectrogram  at  the  top 
was  Blade  with  a  grating,  while  that  at  the  bottom  was  made  with  a  three  prism  spectro- 
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TABLE  I 

VICTORIA  plate,  micrgphotometer,  and  INTENSITOMETER 
DATA  FOR  0  URSAE  MAJORIS 


Victoria 
Plate  . 

.  Number 

Microphotometer  & 

.  Intensitometer 
Tracing  Number 

Spectrograph 

Wavelength  Range 
(Angstroms) 

55224 

1808 

Grating*- 2nd 

Order 

3900-4080 

55053 

1828 

Grating*- 2nd 

Order 

3990-4085 

38157 

1827 

Grating  -3rd 

Order 

4000-4210 

50092  • 

1797 

Prism 

4120-4250 

55225 

1793 

Grating*- 2nd 

Order 

4180-4315 

55054 

1811 

Grating*- 2nd 

Order 

4210-4340 

50092 

1796 

Prism 

4250-4580 

34599 

1807 

Grating  -3rd 

Order 

4275-4455 

55225 

1792 

Grating*- 2nd 

Order 

4300-4480 

37112 

1914 

Grating*- 2nd 

Order 

4310-4900 

37075 

1915 

Grating*- 2nd 

Order 

4310-4900 

55054 

1810 

Grating*- 2nd 

Order 

4335-4495 

31421 

1799 

Prism 

4600-5020 

38133 

1816 

Grating-  -2;nd 

Order 

4700-4830 

38133 

1815 

Grating  -2nd 

Order 

4825-5180 

31421 

1798 

Prism 

5015-5710 

37111 

1812 

Grating  -3rd 

Order 

5130-5450 

37074 

1800 

Grating  -2nd 

Order 

5160-5450 

34800 

1814 

Grating  -2nd 

Order 

5165-5510 

36796 

1818 

Grating  -2nd 

Order 

5600^-6010 

34799 

1817 

Grating  -2nd 

Order 

6000-6325 

36795 

1819 

Grating  -2nd 

Order 

6200-6700 

*The  asterisk  refers  to  grating  spectra  made  with  the  Bausch  and 
Lomb  grating  No.  496,  or  169  in  the  cases  of  plates  37112  and  37075. 
All  other  grating  spectra  were  made  with  the  Wood  grating.  Prism 
spectrograms  were  made  with  the  three  prism  spectrograph. 
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Location  of  the  Continuum 

For  the  most  part,  the  continuum  was  drawn  directly  on  the  intensi- 
tometer  tracings  as  the  average  of  the  galvanometer  deflections  due  to 
plate  grain  in  the  regions  between  the  lines.  In  a  few  cases,  however, 
the  position  of  the  continuum  was  more  apparent  on  the  microphotometer 
traqing  and  hence  was  first  drawn  there  and  then  transferred  to  the 
intensitometer  tracing. 

Approximation  of  the  Profiles 

Theory  indicates  that  the  shape  of  the  wings  exhibited  by  the 
stronger  lines  should  be  inversely  proportional  to  the  square  of  vy, 
the  distanea  from  the  line  center  in  units  of  the  Doppler  width.  An 
attempt  was  made  to  take  this  into  account  when  drawing  the  profiles  of 
the  stronger  lines,  The  shapes  of  the  wings  were  readily  apparent  for 
very  strong  lines,  expecially  when  there  were  no  nearby  lines  to  produce 
blending  effects  in  the  wings.  When  blending  was  serious  the  wings 
were  roughly  approximated . 

Identification  of  the  Lines 

The  lines  were  identified  by  referring  to  the  tables  published  by 
Swensson  (1946)  in  his  paper  "The  Spectrum  of  Procyon."  Identifications 
could  be  made  in  this  way  because  the  spectrum  of  Procyon  (spectral 
type  F5,  luminosity  class  IV)  is  quite  similar  to  that  of  6  Ursae 
Majoris,  Also,  the  0  Ursae  Majoris  tracings  were  compared  with  similar 
ones  for  Procyon  which  were  used  by  Schroeder  (1958)  in  his  study  of 
that  star.  The  work  by  Charlotte  E.  Moore  (1945),  A  Multiplet  Table 
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of  Astrophvsical  Interest.  Revised  Edition,  was  also  consulted  for 
purposes  of  line  identification. 

Selection  of  Lines 

The  selection  of  lines  to  be  used  in  this  study  was  based  upon  the 
following  considerations;  effects  of  blending,  the  availability  of 
f-’values,  and  the  number  of  lines  available  for  a  particular  element. 
These  are  discussed  in  detail  in  the  paragraphs  to  follow. 

Effects  of  Blending 

Equivalent  widths  were  measured  only  for  lines  which  were  believed 
to  be  substantially  free  of  blending  effects  in  the  central  portions  of 
the  profiles.  These  effects  were  tolerated  when  present  in  the  wings 
of  the  strong  lines,  since  this  was  almost  always  the  case.  Decisions 
as  to  whether  or  not  blending  was  serious  enough  to  preclude  the  use 
of  a  line  were  based  upon  Swensson's  line  identifications  for  Procyon. 
A,ny  line  listed  by  Swensson  as  having  more  than  one  contributor  was 
discarded.  Also,  the  selected  lines  were  restricted  to  those  with 
wavelengths  greater  than  4000  angstroms.  To  the  violet  of  X.4000  the 
spectrum  is  very  complex,  and  many  lines  in  this  region  are  on  the 
broad  wings  of  the  H  and  K  lines  of  ionized  calcium. 

Availability  of  f- Values 

An  extensive  bibliography  on  transition  probabilities  has  been 
published  by  G!l  ennon  and  Wiese  (1962).  This  paper  proved  to  be  of 
valuable  assistance  in  gathering  information  concerning  f-values. 


31 


Wright  (1948)  has  pointed  out  that  oscillator  strengths  measured 
in  the  laboratory  are  more  accurate  than  those  calculated  theoretically. 
Corliss  and  Bowman  (1962)  have  published  laboratory  absolute  f-values 
for  25,000  lines  between  2000  and  9000  angstroms  for  112  spectra  of  70 
elements.  Their  work  was  done  at  the  National  Bureau  of  Standards,  and 
henceforth  these  f-values  will  be  referred  to  as  "NBS  f-values."  The 
existence  of  this  extensive  work  should  greatly  reduce  the  use  of 
thepretically  calculated  f-values  in  studies  of  the  type  presented  in 
this  paper.  However,  some  limitations  are  present  in  these  listings. 

Due  to  the  relatively  low  temperature  (about  5100°K)  of  the  copper  arc 
used  to  excite  the  lines,  the  weaker  lines  of  several  of  the  spectra, 

Ti  II  for  example,  were  not  excited  sufficiently  for  measurement. 

Henpe,  no  f-values  are  given  for  these  lines.  The  accuracy  of  the  NBS 
f-values  is  admittedly  not  as  great  as  for  f-values  measured  by  more 
refined  and  laborious  techniques,  but  the  NBS  measurements  are  almost 
certainly  more  reliable  than  theoretical  f-values.  The  principal 
usefuiness  of  the  NBS  tables  lies  in  the  large  number  of  values  made 
available,  the  wide  range  of  wavelength  covered,  and  the  fact  that  the 
scale  is  absolute. 

Tables  of  laboratory  oscillator  strengths  for  Fe  I  and  Ti  I  have 
been  published  by  the  Kings  (1938),  for  another  group  of  Fe  I  lines 
by  Carter  (1949),  and  for  Cr  I  by  Hill  and  King  (1951).  All  of  these 
measurements  are  considered  to  be  quite  accurate.  These  f-values  are 
on  relative  scales. 

All  of  the  above  mentioned  works  were  used  as  sources  of  f-values, 
although  the  NBS  tables  were  the  major  source.  Most  of  the  lines  for 
which  f-values  are  given  by  the  Kings,  Carter,  and  Hill  are  also  included 
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in  the  NBS  tabulations;  but  due  to  the  more  accurate  nature  of  the 
measures  of  the  Kings,  Carter,  and  Hill,  it  was  considered  advisable 
to  make  independent  analyses  of  the  curves  of  growth  for  Fe  I,  Ti  I, 
and  Cr  I  using  these  oscillator  strengths.  The  results  obtained  for 
the  different  sets  of  f-values  could  then  be  compared.  Also,  the  Kings' 
and  Carter's  f-values  for  Fe  I  were  not  combined,  but  were  used  to 
derive  independent  results. 

As  mentioned  above,  there  was  found  to  be  a  shortage  in  the  NBS 
tables  of  f-^values  for  the  weaker  lines  of  Ti  II.  This  was  the  only 
ionized  element  considered  suitable  for  investigation  here,  and  thus 
offered  the  only  opportunity  to  compare  abundances  of  the  neutral  and 
ionized  atoms.  To  improve  the  accuracy  of  the  results  obtained  for 
Ti  II,  it  was  decided  to  supplement  the  available  data  by  including  a 
number  of  weaker  lines  for  which  Wright  (1948)  has  tabulated  solar 
log  X^-values.  These  can  be  converted  to  f-values  in  the  manner  der 
scribed  in  Wright's  paper,  i.e.,  by  means  of  the  relation 

log  gf\  =  log  +  5040  Xg/T,  (22) 

where  Xg  is  ti*®  excitation  potential  in  electron  volts  of  the  lower 
level  of  the  transition  in  question,  and  T  =  4700°K  is  a  mean  value, 
suggested  by  Wright,  of  the  excitation  temperature  for  the  sun.  The 
application  of  this  equation  puts  the  derived  f-values  on  a  relative 
scale.  Values  were  calculated  for  ten  lines  of  Ti  II  which  are  common 
to  Wright's  and  the  NBS  tabulations,  and  the  average  value  of  the 
differences  between  these  and  the  corresponding  NBS  quantities  was 
used  to  convert  to  an  absolute  scale,  that  of  the  NBS.  Figure  5 
compares  the  derived  oscillator  strengths  (on  the  absolute  scale)  with 
the  NBS  values  for  these  common  lines. 


The  reliability  of  the  results  obtained  by  the  method  of  curve 
of  growth  analysis  employed  in  this  study  depends  strongly  upon  the 
number  of  lines  available  for  the  element  being  investigated.  Hence, 
it  was  felt  that  results  derived  from  only  a  few  lines  would  be  of 
doubtful  value.  Of  the  elements  studied,  Co  I  was  represented  by 
twelve  lines.  Elements  with  fewer  lines  available  were  not  considered 
foy  analysis. 

Table  II  lists  the  eight  atoms  and  ions  studied,  the  sources 
providing  the  f- values,  and  the  number  of  lines  available  for  each  set 


of  f-f values. 
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TABLE  II 

LIST  OF  ATOMS  AND  IONS  STUDIED  IN  THIS  INVESTIGATION 


Atom 
or  Ion 

Source 
of  f-values 

No .  of 

Lines 

Atom 
or  Ion 

Source 
of  f-values 

No .  of 
..Lines 

Ga  I 

NBS 

21 

Mn  I 

NBS 

22 

Ti  I 

NBS 

55 

Fe  I 

NBS 

115 

Ti  X 

King 

17 

Fe  I 

King 

30 

Ti  11 

NBS  &  Wright 

27 

Fe  I 

Carter 

29 

Cr  I 

NBS 

65 

Co  I 

NBS 

12 

Cr  I 

Hill  &  King 

36 

Ni  I 

NBS 

32 

Determination  of  Equivalent  Widths 

The  determination  of  equivalent  widths  first  involves  measurement 
of  the  area  enclosed  by  the  line  profile  and  the  continuum.  These  areas 
were  obtained  using  an  Ott  rolling  disc  planimeter.  Each  line  was 
measured  twice,  each  measurement  consisting  of  two  circuits  of  the 
profile  with  the  planimeter.  The  continuum  height  was  taken  as  the 
average  of  the  values  adjacent  to  the  line.  For  each  tracing  the  dis¬ 
persion  was  taken  at  various  points  along  the  spectrum  and  a  straight 
line  fitted  to  a  plot  of  dispersion  against  wavelength.  For  the  calcu¬ 
lation  of  equivalent  widths  the  dispersion  was  then  read  from  these 
plots,  W  was  calculated  in  milliangstroms .  Due  to  overlapping  and 
duplication  of  the  spectrograms,  some  lines  were  represented  by  as 
many  as  seven  profiles,  a  few  by  only  one.  The  equivalent  width,  in 
all  cases,  was  taken  as  the  average  with  all  measures  assigned  equal 
weight,  For  a  representative  group  of  lines.  Figure  6  compares  the 
equivalent  widths  measured  on  different  tracings  with  the  corresponding 
average  values.  Every  tracing  used  in  this  study  is  represented  in 
this  plot  by  at  least  one  measurement  of  W. 


log  W  (mA) 
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log  W  (mA) 
avg 

Figure  6.  Comparison  of  Equivalent  Widths  Measured  on  Different  Tracings 
•With  the  Average  Values  for  a  Representative  Group  of  Lines.  The  scale 
is  logarithmic  with  the  average  values  as  abscissae  and  the  values 
obtained  from  different  tracings  as  ordinates.  A  set  of  points  forming 
a  vertical  array  represents  the  values  of  W  measured  for  a  particular 
line.  The  45-degree  line  gives  the  position  of  the  average. 


36 


Thackeray's  (1936)  relation  was  applied  to  increase  the  measured 
equivalent  widths  of  the  lines  located  in  the  wings  of  the  Balmer  lines, 
This  relation  is 


W  = 


(23) 


w 

where  is  fhe  measured  equivalent  width  referred  to  the  wing  as  the 
continuum,  r^  is  the  ratio  of  the  intensity  of  the  wing  to  the  intensity 
pf  the  true  continuum  at  the  center  of  the  line  being  studied,  and  W 
fs  the  equivalent  width  the  line  would  have  if  no  blending  occurred. 

This  relation  is  not  exact  but  should  be  sufficiently  accurate  to 
permit  the  use  of  such  lines. 


Table  of  Equivalent  Width  Measures 


Equivalent  widths  obtained  in  this  study  are  listed  in  Table  III 
along  with  other  quantities  used  in  the  analysis.  A  section  is  included 
for  each  element  studied,  and  these  are  arranged  according  to  increasing 
atomic  weight.  Within  each  section  the  lines  are  tabulated  in  the  order 
of  increasing  wavelength. 

Column  1  lists  the  wavelength  in  angstroms  as  given  by  Miss  Moore 
(1945)  in  the  Revised  Multiplet  Table  (RMT)  . 

Column  2  gives  the  RMT  multiplet  number. 

Column  3  lists  the  excitation  potential  of  the  lower  level  of  the 
transition,  in  electron  volts. 

Column  4  gives  log  gfX,  where  the  wavelength  is  in  centimeters. 


This  column  is  divided  into  two  sub- columns.  The  sub- column  on  the 
left,  headed  by  "NBS,"  contains  values  of  log  gfX  taken  from  the  NBS 
tables.'  The  sub- column  on  the  right,  headed  by  "Other,"  contains 
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values  of  log  gfA.  from  other  sources;  for  Tf- I  the  values  in  this 
sub-' column  are  taken  from  the  Kings'  tables,  for  Ti  II  the  values  are 
from  Wright,  for  Cr  I  the  values  are  from  Hill  and  King,  and  for  Fe  I 
the  values  are  from  the  Kings'  tables  and  Carter's  tables.  Carter's 
values  are  followed  by  the  letter  "C"  so  that  values  from  the  two 
sources  may  be  distinguished. 

Column  5  lists  the  number  of  profiles  measured  for  the  line. 

Column  6  gives  the  equivalent  width  in  milliangstroms .  Equivalent 
widths  of  lines  found  in  the  wings  of  the  Balmer  lines  have  been  correct¬ 
ed  according  to  Thackeray's  relation  and  are  indicated  by  an  asterisk. 
Data  pertinent  to  these  lines  are  listed  in  Table  IV. 

Column  7  gives  log  W/\,  where  W  and  X  have  the  same  units. 

Column  8  lists  the  value  of  the  limb  darkening  correction,  Alog  T|^, 

for  the  Milne- Eddington  pure  scattering  model.  Such  corrections  can 

be  obtained  by  interpolating  in  the  published  curve  of  growth  tables. 

The  addition  of  Alog  T]^  to  the  abscissa  of  each  line  accounts  for 

differences  in  limb  darkening  (which  depends  upon  the  wavelength)  and 

makes  possible  the  fitting  of  all  lines  to  a  theoretical  curve  of  growth 

0,1 

calculated  for  a  single  value  of  B  /B  ,  here  taken  to  be  2/3. 

Alog  T|^  is  a  function  of  B°/B^.  Hence,  to  determine  Alog  T)^  for 
a  given  line,  it  is  necessary  to  know  the  value  of  B°/B^  corresponding 
to  the  wavelength  of  the  line.  It  was  assumed  that  the  values  of 
B°/B^  for  0  Ursae  Majoris  are  approximately  the  same  as  for  the  sun, 
where  the  source  of  continuous  opacity  is  the  negative  hydrogen  ion. 

This  seems  to  be  a  valid  assumption,  inasmuch  as  the  temperature  of 
9  Ursae  Majoris  is  only  slightly  higher  than  that  of  the  sun  and  is 
certainly  much  too  low  for  neutral  hydrogen  and  helium  to  play  a 
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DATA  PERTINENT  TO  LINES  LOCATED  ON  WINGS  OF  HYDROGEN  LINES 


Element 


Wing 

Location 


4092.386 

Co 

I 

H6 

4107.492 

Fe 

I 

H6 

4110.532 

Co 

I 

H6 

4337.049 

Fe 

I 

Hy 

4337.566 

Cr 

I 

Hy 

4337.916 

Ti 

II 

Hy 
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Cr 

I 

Hy 
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Cr 

I 

Hy 
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Ti 

II 

Hy 
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Cr 

I 

Hy 
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Cr 

I 

Hy 
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Fe 

I 

HP 

4866.267 

Ni 

I 

HP 

4871.323 

Fe 

I 

HP 

4872.144 

Fe 

I 

HP 

4873.437 

Ni 

I 
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Corrected 
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112 

71 

85 

34 

38 

27 

158 

71 

93 

17 
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30 

54 

13 

25 

84 

104 

81 

97 

26 

29 

87 

145 

54 

66 

74 

155 

72 

191 

48 
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prominent  role  in  the  continuous  absorption.  Therefore,  the  values  of 
B°/B^  used  in  this  paper  are  taken  from  the  solar  values  observed  by 
Houtgast  (1942) .  The  values  were  read  from  the  plot  in  Figure  7  which 
was  constructed  from  Houtgast 's  data. 

Alog  T|^  varies  with  the  intensity  of  the  line,  i.e..,  with  log  W/b, 
as  well  as  with  B°/B^.  In  order  to  obtain  approximate  values  of  log  W/b 
for  each  line,  log  c/v  was  taken  to  be  approximately  equal  to  5.00; 
values  in  this  neighborhood  had  been  obtained  from  the  applications  of 
the  other  models.  Then,  using  the  measured  values  of  log  W/\, 


log  W/b  =  log  W/\  +  5.00. 


(24) 


The  limb  darkening  corrections  are  slightly  affected  by  the  value 


of  the  damping  parameter,  a.  These  variations  were  considered  to  be 


negligibly  small,  however,  and  the  values  of  Alog  7]^  given  in  Column  8 
are  for  log  a  =  -1.4. 
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Figure  7,  B°/B  Versus  Wavelength,  The  data  used  to  construct  the 
curve  are  the  solar  values  observed  by  Houtgast  (1942) .  The  circled 
points  are  the  values  given  by  Houtgast.  The  points  were  connected 
by  a  smooth  curve  for  purposes  of  interpolation, 
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For  the  S'-S  pure  scattering  model,  limb  darkening  corrections  are 
not  large  enough  to  be  of  any  conpequence  and  have  not  been  applied  in 
this  study. 

Limb  darkening  effects  are  taken  into  account  in  a  different  way 
by  Jlunger  in  his  pure  absorption  curves.  Here,  the  limiting  depth,  R^, 
which  depends  upon  appears  in  the  ordinate  of  the  theoretical 

curves.  The  observed  curves  of  growth  in  the  case  of  pure  absorption 
are  plotted  with 

log  W/2R  X  =  log  W/X  -  log  2R  (25) 

c  c 

as  ordinate.  The  values  of  log  2R^  were  read  for  each  line  from 
Table  V,  which  gives  log  2R^  as  a  function  of  B°/B^.  In  order  to  form 
this  table,  R^  was  calculated  from  Equation  (5)  for  each  value  of  B°/B^ 
listed . 


Comparison  of  Equivalent  Width  Measurements 
With  Those  of  Greenstein 

In  Figure  8,  a  comparison  is  made  of  the  equivalent  widths  of  the 
9  Ursae  Majoris  lines  measured  in  this  paper  and  identical  lines  measured 
by  Greenstein  (1948).  Inspection  of  the  figure  reveals  good  agreement 
for  values  of  log  W/X  around  -4.50.  For  stronger  lines,  our  values  are 
larger  than  those  of  Greenstein,  and  the  reverse  is  true  for  the  weaker 
lines.  These  variances  may  be  due  to  the  different  methods  used  to 
obtain  equivalent  widths.  In  Greenstein' s  work  the  prpfiles  were  drawn 
on  the  microphotometer  tracings,  line  depths  measured  every  millimeter, 
converted  into  absorptions,  and  numerically  integrated  to  obtain  the 
equivalent  widths.  In  this  paper  the  profiles  were  drawn  on  the  intensi- 
tometer  tracings  rather  than  the  microphotometer  tracings.  The  profiles 
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were  then  planimetered  to  obtain  equivalent  widths.  Consequently, 
differences  might  be  expected. 


TABLE  V 

LOG  2R  AS  A  FUNCTION  OF  B°/B^ 
c 


log  2R 

c 

b°/b‘ 

log  2R^ 

B°/B^ 

log  2R^ 

0.20 

+0 . 19 

0.43 

+0.08 

0.66 

0.00 

0,21 

+0.18 

0.44 

+0.08 

0.67 

0.00 

0,22 

+0.18 

0.45 

+0.08 

0.68 

0.00 

0.23 

+0.17 

0.46 

+0.07 

0.69 

-0.01 

0.24 

+0.17 

0.47 

+0.07 

0.70 

-0.01 

0.25 

+0,16 

0.48 

+0.07 

0.71 

-0.01 

0.26 

+0.16 

0.49 

+0.06 

0.72 

-0.02 

0.27 

+0.15 

0.50 

+0.06 

0.73 

-0.02 

0,28 

+0.15 

0.51 

+0,05 

0.74 

-0.02 

0.29 

+0.14 

0.52 

+0.05 

0.75 

-0.03 

0.30 

+0.14 

0.53 

+0.05 

0.76 

-0.03 

0.31 

+0.14 

0.54 

+0.04 

0.77 

-0.03 

0.32 

+0.13 

0.55 

+0.04 

0.78 

-0.04 

0.33 

+0.13 

0.56 

+0.04 

0.79 

-0.04 

0.34 

+0.12 

0.57 

+0.03 

0.80 

-0.04 

0.35 

+0.12 

0.58 

+0.03 

0.81 

-0.04 

0.36 

+0.11 

0.59 

+0.03 

0.82 

-0.05 

0.37 

+0.11 

0.60 

+0.02 

0.83 

-0.05 

0.38 

+0.10 

0.61 

+0.02 

0.84 

^0.05 

0,39 

+0.10 

0.62 

+0,02 

0.85 
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0.40 
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0.63 

+0.01 
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-0.06 

0.41 
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0.88 
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Figure  8.  Comparison  of  Equivalent  Widths  for  Lines  Common  to  This  Paper 
and  Greenstein. 


CHAPTER  IV 


.  RESULTS 

Discussion  of  the  Procedures  Used  in  the  Reduction  of  the  Data 

The  general  procedures  employed  in  this  study  to  determine  log  c/v, 
log  a,  the  excitation  temperature,  and  the  abundance  from  the  curve  of 
growth  have  been  discussed  in  Chapter  II,  The  details  involved  in  the 
application  of  these  procedures  are  considered  in  the  following  para¬ 
graphs,  For  convenience  in  the  discussion,  the  designations  "Fe  I-NBS," 
"Fe  I-King,"  and  "Fe  I-Carter"  are  used  to  differentiate  between  the 
three  different  sets  of  data  for  Fe  I  corresponding  to  the  three  sources 
of  f- values.  Similarly,  data  pertaining  to  Ti  I  and  Cr  I  are  referred 
to  as  "Ti  I-NBS,"  "Ti  I-King,"  "Cr  I-NBS,"  and  "Cr  I^Hill"  data. 

Determination  of  Log  c/v 

To  provide  a  means  for  Comparing  the  observed  curve  of  growth  for 
the  lines  of  a  given  element  in  a  particular  stage  of  ionization  with 
the  theoretical  curves,  the  values  of  log  W/X  (or  log  W/2R^X)  for  lines 
which  all  arise  from  a  given  atomic  term  were  plotted  against  the  corres' 
ponding  values  of  log  gf\  (or  log  gfX  +  Alog  7]^) .  Such  plots  represent 
portions  of  the  observed  curve  of  growth.  Each  point  on  these  plots 
was  weighted  roughly  according  to  the  number  of  profiles  measured  in 
obtaining  the  equivalent  width  of  the  line-- the  more  profiles  measured, 
the  more  reliable  the  point.  The  sizes  of  the  points  were  made  to 
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indicate  the  weights,  and  the  larger  (more  reliable)  points  were  given 
greater  influence  in  determining  the  fit  of  the  plot  on  the  theoretical 
curve  of  growth. 

It  is  necessary  to  have  a  plot  which  defines  the  knee  of  the  observed 
curve  of  growfh  in  order  to  determine  log  c/v.  The  vertical  shift  re- 

I 

quired  to  fit  such  a  plot  on  the  theoretical  curve  gives  log  c/v,  as  can 
be  seen  from  Equation  (16)  of  Chapter  II.  In  this  manner,  values  were 
obtained  from  the  observed  curves  for  Ti  I-NBS,  Ti  I-King,  Ti  II,  Cr  I- 
NBS,  Cr  I“Hill,  Fe  I-NBS,  and  Fe  I-King.  In  each  case,  four  values  of 
log  c/v  were  determined  corresponding  to  the  four  different  atmospheric 
models.  No  plots  suitable  for  the  determination  of  log  c/v  were  avail¬ 
able  for  any  of  the  other  sets  of  lines  investigated, 

Because  of  the  large  number  of  lines  available  and  the  relatively 
small  scatter  of  the  points  defining  the  observed  curve  of  growth,  the 
values  of  log  c/v  obtained  for  Fe  I-NBS  were  regarded  as  more  reliable 
than  any  of  the  other  determinations.  Hence,  the  Fe  I-NBS  values  of 
log  c/v  were  adopted  for  use  with  the  corresponding  atmospheric  models 
in  the  analysis  of  the  lines  of  Ca  I,  Mn  I,  Co  I,  and  Ni  I.  For  Fe  I- 
Carter,  the  one  remaining  set  of  lines,  the  values  of  log  c/v  found  for 
Fe  I-King  were  employed.  Since  the  Kings'  f-values  were  used  by 
Carter  to  establish  a  preliminary  relative  scale  for  his  measurements, 
these  two  sets  of  f-values  are  closely  related.  For  this  reason,  the 
Fe  I-King  values  of  Iqg  c/v  were  felt  to  be  the  most  appropriate  choices 
for  use  with  the  Fe  I-Carter  lines. 

Concerning  the  validity  of  using  the  Fe  I-NBS  values  of  log  c/v 
in  the  study  of  the  lines  of  Ca  I,  Mn  I,  Co  I,  and  Ni  I,  it  should  be 
noted  that  the  velocity,  v,  depends  upon  the  mass  of  the  atom  in  question. 
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Howevey,  this  dependence  should  be  negligibly  small,  as  indicated  by 
the  following  discussion.  Equation  (2)  of  Chapter  II  can  be  written  as 

■  \l  ''turb  ■  «6) 

where  M  is  the  mass  of  the  atom  and  T  is  the  kinetic  temperature. 

Before  this  equation  can  be  used  to  compare  values  of  v  for  different 
atoms,  it  is  necessary  to  specify  a  kinetic  temperature  and  a  turbulent 
velocity.  It  was  assumed  that  the  kinetic  temperature  is  approximately 
equal  to  the  effective  temperature,  i.e.,  T  =  6210°K  as  listed  by  Keenan 
and  Morgan  (1951)  for  stars  of  class  F6  IV.  Then,  using  the  mass  of 
the  iron  atom  along  with  the  Fe  I-NB3  values  of  log  c/v,  turbulent 
velocities  were  calculated  from  Equation  (26) .  The  four  values  obtained, 
corresponding  to  the  four  different  atmospheric  models,  were  found  to 
range  between  2.16  x  10^  and  2.67  x  10^  cm/sec.  Using  these  estimated 
values  of  Equation  (26)  was  employed  to  calculate  values  of  v  for 

calcium,  which  has  an  atomic  mass  differing  from  that  of  iron  by  a 
greater  amount  than  any  of  the  other  elements  investigated,  It  was 
found  that  the  resulting  values  of  log  c/v  were  substantially  the  same 
as  those  for  Fe  I-NBS.  For  example,  the  largest  difference  between 
the  Fe  I-NBS  and  calcium  values  for  the  same  atmospheric  model  amounts 
to  0,02,  which  is  negligible  for  the  purposes  of  this  investigation. 
Because  of  these  results,  it  was  felt  that  the  values  of  log  c/v  deter¬ 
mined  for  Fe  I-NBS  were  the  best  approximations  available  when  this 

•. 

quantity  could  not  be  determined  directly. 

Determination  of  Log  a 

*  If  strong  lines  defining  the  damping  portion  of  the  curve  of  growth 
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are  available,  log  a  can  be  determined.  For  Fe  I-NBS,  Fe  I-King,  and 
Ti  II,  single  plots  for  a  given  term  were  available  with  enough  lines 
of  a  sufficient  range  of  intensity  to  adequately  define  the  transition 
portion  (the  section  between  the  knee  and  the  beginning  of  the  damping 
portion)  as  well  as  the  damping  portion.  Since  the  transition  portion 
of  the  curve  of  growth  is  not  perfectly  horizontal,  the  onset  of  damp¬ 
ing  occurs  at  greater  intensities  as  the  damping  parameter  decreases. 
Hence,  fitting  the  above  mentioned  plots  on  the  theoretical  curves  for 
different  values  of  the  damping  parameter  requires  larger  vertical 
shifts  for  smaller  values  of  log  a.  However,  the  vertical  shifts  were 
determined  by  other  plots  that  fit  around  the  knee.  In  each  case,  using 
the  appropriate  vertical  shift,  the  plot  defining  the  transition  and 
damping’’ portions  was  superposed  upon  theoretical  curves  of  growth  for 
different  values  of  log  a.  The  value  belonging  to  the  curve  giving 
the  best  fit  was  adopted  as  characteristic  of  the  element  and  atmospheric 
model  under  investigation. 

For  Ca  I,  Cr  I-NBS,  Cr  I-Hill,  and  Mn  I  only  a  few  lines  were 
strong  enough  to  fit  on  the  damping  portion.  Since  it  was  not  possible 
to  locate  the  position  of  the  change  from  transition  to  damping  with 
the  plots  available  here,  the  procedure  discussed  in  the  above  para" 
graph  could  not  be  used.  For  each  of  these  cases  the  value  of  log  a 
was  taken  to  be  that  yhich  resulted  in  the  best  straight  line  when 
Alog  X  was  plotted  against  .  This  straight  line  relation  is  discussed 
further  in  the  next  section  which  is  concerned  with  the  determination  of 

"it 

Since  the  excitation  potentials  of  the  lower  levels  of  the  transi¬ 
tions  may  be  slightly  different  for  different  lines  arising  from  a 
particular  atomic  term,  an  average  value,  was  used  for  each  term. 
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excitation  temperatures.  A  set  of  values  of  Alog  X  was  obtained  for 
each  value  of  log  a,  and  Alog  X  plotted  against  .  The  value  of  log  a 
giving  the  best  straight  line  was  determined  by  comparing  these  plots. 
This  rather  indirect  method  for  determining  log  a  should  not  be  consid¬ 
ered  as  reliable  as  the  procedure  described  in  the  preceding  paragraph. 

Because  of  the  absence  of  strong  lines,  it  was  not  possible  to 
determine  values  of  log  a  from  the  lines  of  Ti  I-NBS,  Ti  I-King,  Fe  I- 
Carter,  Co  I,  and  Ni  I.  For  a  particular  atmospheric  model,  all  the 
theoretical  curves  of  growth  for  various  values  of  log  a  are  the  same 
below  the  knee.  Hence,  the  value  of  log  a  is  of  no  significance  when 
only  weak  lines  are  considered. 

Determination  of  the  Excitation  Temperature 

Equation  (18)  of  Chapter  II  may  be  written  as 

Alog  X  =  constant  -  (27) 

where  v  is  the  average  of  the  excitation  potentials  of  all  the  lines 
arising  from  a  given  term.  It  is  apparent  from  this  equation  that  a 
straight  line  relation  should  exist  between  Alog  X  and  with  a  slope 
equal  to  -9.  For  each  atomic  term  Alog  X  was  found  from  the  hori¬ 
zontal  shift  required  to  fit  the  corresponding  plot  on  the  theoretical 
curve  of  growth.  The  same  vertical  shift  was  used  for  fitting  to  a 
particular  curve  all  plots  belonging  to  the  same  element  in  a  given 
stage  of  ionization,  although  this  vertical  shift  usually  changed  in 
going  from  one  atmospheric  model  to  another.  The  method  of  least  squares 
was  employed  to  obtain  the  slope,  -0,  of  the  straight  line  which  best 
fit  the  plot  of  Alog  X  versus  x^.  The  excitation  temperature  was  then 
calculated  from  T  ■  5040/9. 
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For  the  purposes  of  the  least  squares  calculations  each  value  of 
Alog  X  was  assigned  a  weight  indicative  of  its  apparent  reliability.  A 
weight  was  assigned  by  considering  the  number  of  points  on  the  plot 
used  to  determine  Alog  X,  the  scatter  of  the  points,  the  extent  of  the 
curve  of  growth  defined  by  the  plot,  and  the  fit  of  the  plot  on  the 
theoretical  curve  of  growth.  Some  terms  were  represented  by  only  one 
point.  In  such  cases  the  values  of  Alog  X  obtained  were  given  zero 
weight  when  there  were  a  sufficient  number  of  other  values  to  adequately 
determine  the  straight  line  relating  Alog  X  and  .  Also,  when  a  value 
of  Alog  X  was  clearly  far  off  the  preliminary  straight  line,  it  was 
given  zero  weight. 

Determination  of  the  Abundance 

In  Chapter  II  the  equation  for  the  abundance  in  the  case  of  S-S 
pure  scattering  is  presented  (Equation  (20)  ) ,  and  the  means  of  modifying 
this  equation  so  that  it  will  apply  to  any  of  the  other  models  is  indi-- 
cated.  Written  out  in  full,  the  equations  giving  the  abundances  fpr 
the  four  models  are  as  follows; 

M-E  pupe  scattering  and  pure  absorption^- 

log  W/pH  ^  shift  -  log  C  +  log  u  log  c/v,  (28) 

S-S  pure  scattering-- 

log  NH  *  shift  -  log  C  +  log  u  -  log  c/v,  (29) 

and  S-S  pure  absorption-- 

log  NH  ■  shift  "  log  C  +  log  u  ■?  log  c/v  +  log  3/2.  (30) 

When  relative  f-values  were  used,  as  was  the  case  for  Ti  I-King, 

Cr  I-Hilli  Fe  I-King,  and  Fe  I-Carter,  it  was  necessary  to  subtract 
another  quantity,  Alog  gf,  from  the  right-hand  sides  of  Equations  (28), 
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(29),  and  (30).  Alog  gf  is  defined  by  the  following; 

ilog  g£  -  log  -  log  .  (31) 

The  use  of  relative  f-values  causes  the  shifts  obtained  to  be  too  large 
by  the  amount  Alog  gf;  hence  the  subtraction  of  this  quantity.  If  one 
or  more  absolute  f-values  are  known  for  lines  included  in  the  tables  of 
relative  fi-values,  which  are  on  a  uniform  scale,  Alog  gf  can  be  cal¬ 
culated  from  Equation  (31),  The  values  used  for  Ti  I-Kfng,  Cr  I-Hill, 

Fe  iTKing,  and  Fe  I-Carter  are  given  later  in  this  chapter  when  the 
results  obtained  from  these  sets  of  lines  are  presented. 

As  pointed  out  in  Chapter  II,  shift  =  Alog  X  +  From  this 

equation,  values  of  the  shift  were  determined  for  each  atomic  term,  and 
weighted  average  values  were  calculated  using  the  same  weights  as 
employed  in  the  least  squares  determinations  of  the  excitation  temper¬ 
atures.  The  weighted  average  values  were  used  in  the  abundance  compu¬ 
tations.  The  constant  C  is  given  by  Equation  (9)  of  Chapter  II.  The 
values  of  log  u  were  obtained  by  interpolating,  when  necessary,  in  the 
tables  published  by  Aller  (1960),  where  log  u  is  given  as  a  function  of 
the  excitation  temperature. 

Values  have  been  determined  in  this  study  for  the  abundances  in 
the  atmosphere  of  9  Ursae  Majoris  of  the  neutral  atoms  of  seven  elements 
and  the  singly  ionized  atoms  of  one  element,  titanium.  Titanium  was 
the  only  element  studied  where  results  were  obtained  for  both  the  neutral 
and  singly  ionized  atoms.  Since  second  and  higher  ionizations  are 
negligible,  the  total  abundance  of  titanium  can  be  obtained  by  adding 
the  abundances  of  the  neutral  and  singly  ionized  atoms.  It  was  found 
that  the  contribution  of  the  neutral  atoms  to  the  total  abundance  is 
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completely  negligible  as  far  as  the  accuracy  attained  in  this  study  is 
concerned.  Hence,  the  abundances  obtained  for  singly  ionized  titanium 
also  represent  the  total  abundances  for  titanium. 

For  the  other  six  elements  studied,  where  only  the  abundances  of 
the  neutral  atoms  were  determined,  the  Saha  equation  provides  a  means 
for  estimating  the  contributions  of  the  ionized  atoms  to  the  total 
abundances.  The  validity  of  the  results  obtained  is  affected  by  depar¬ 
tures  ftom  idealized  conditions,  Although  such  departures  (e.g.,  non- 
uniform  distribution  of  the  atoms  and  ions  throughout  the  stellar 
atmosphere  and  deviations  from  thermodynamic  equilibrium)  are  known  to 
exist,  the  Saha  equation,  being  the  best  means  available,  was  applied 
to  estimate  the  total  abundances. 

The  Saha  equation  is  presented  in  Chapter  II,  and  for  convenience 
it  is  restated  below: 

(NH)  2uj^ 

log  =  -log  ^  Xi  +  2.5  log  T  .  0,48  +  log  (32) 

0  o 

As  stated  here  the  equation  applies  to  the  S-S  model.  The  subscripts 
1  and  0  refer  respectively  to  the  singly  ionized  atoms  and  to  the  neutral 
atoms.  For  application  to  the  M-E  model,  (NH)^^  and  (NH)^  are  replaced 
by  (N/phOj^  and  (N/pn)^. 

For  purposes  of  the  calculations,  the  temperature,  T,  was  taken  to 
be  equal  to  the  effective  temperature,  6210°K,  which  is  characteristic  ’ 
of  an  optical  depth  in  the  atmosphere  of  0.6.  This  is  a  representative 
depth  for  line  formation.  It  should  be  noted  that  results  obtained  from 
the  Saha  equation  ate  sensitive  to  the  value  used  for  the  temperature, 
which  leads  to  some  uncertainty  in  the  results.  The  ionization  poten¬ 
tials,  x^»  were  taken  from  the  Revised  Multiolet  Table.  Values  of 
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the  partition  functions,  and  which  depend  upon  the  temperature, 
were  obtained  by  referring  to  the  table?  published  by  Aller  (1960). 

The  one  remaining  quantity  necessary  for  the  calculations  is  the 
electron  pressure,  P^,  The  fact  that  abundances  were  determined  in  this 
investigation  for  both  neutral  and  singly  ionized  atoms  of  titanium 
offered  the  opportunity  to  calculate  values  of  from  the  Saha  equation. 
For  this  purpose,  T  was  taken  to  be  equal  to  6210°K  as  discussed  above, 

^o’  *^1  obtained  from  the  above  mentioned  sources. 

Since  abundances  were  determined  for  four  different  atmospheric  models, 

2 

four  values  of  log  P^,  where  P^  is  in  dynes/cm  ,  were  obtained:  1.36, 
1.30,  1.27,  and  1.32  which  correspond  to  M-E  pure  scattering,  M-E  pure 
absorption,  S-S  pure  scattering,  and  S-S  pure  absorption  respectively. 

The  abundances  of  the  neutral  atoms  determined  for  Ti  I-NBS  were  used 
in  the  calculations.  The  average  of  the  values  obtained  for  log  P 
is  1.31.  This  average  value  was  used  in  subsequent  calculations. 

Since  second  and  higher  ionizations  are  negligible,  the  total 
abundance  in  the  S-S  model  is  given  by 


("^^total  “  ^"^>0  <*^>1 


(NH). 


1  + 


(NH) 

InhT 


11 


(33) 


It  follows  that 


log  (NH)^^^g^  -  log  (NH)^  +  Alog  ?1, 


(34) 


where 

r  ^^hi 

Alog  log  1^1  +  ■(jjjjy'j.  (35) 

Replacing  NH  by  N/pn  makes  the  above  equations  valid  for  the  M^E  model. 
Alog  gives  the  contribution  of  the  ionized  atoms  to  log  (HH)^p^^j^  or 

log  •  For  each  of  the  six  elements  for  which  the  abundances 
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of  the  neutral  atoms  only  were  determined,  the  Saha  equation  was  employed 
to  calculate  a  value  of  (NH)j^/(NH)^  from  which  Alog  7)  was  obtained.  For 
a  given  element  (NH)j^/(NH)^  and  (N/pTi) (N/pn)^  have  the  same  value. 

Hence  Alog  9^  does  not  depend  upon  the  atmospheric  model. 

The  results  of  the  determinations  of  the  values  of  log  c/v  and 
log  a,  the  excitation  temperatures,  and  the  abundances  of  the  atoms  for 
each  set  of  lines  studied  are  presented  in  the  following  paragraphs. 

The  order  of  presentation  is  the  order  in  which  the  sets  of  lines  are 
listed  in  Table  II  of  Chapter  III. 

Results  for  Ca  I 

The  f- values  employed  here  were  those  of  the  NBS .  Observational 
data  were  available  for  twenty-one  lines  belonging  to  five  different 
lower  atomic  terms.  Table  VI  lists  the  terms,  the  values  of  in 
electron  volts,  the  numbers  of  lines  representing  each  term,  the  weights 
used  in  the  least  squares  calculations  of  the  excitation  temperatures, 
and  the  values  of  Alog  X  and  the  shift  obtained  from  the  application 
of  each  of  the  four  sets  of  theoretical  curves  of  growth.  The  values 
of  log  c/v  'listed  in  the  top  row  of  Table  VI  were  adopted  from  those 
obtained  for  Fe  I-iNES.  This  is  indicated  in  the  table  by  enclosing  the 
values  in  parentheses.  The  second  row  gives  the  values  of  log  a,  which 
were  determined  by  the  second  method  discussed  earlier  in  this  chapter 
for  the  determination  of  these  quantities,  i.e.,  from  the  value  giving 
the  best  straight  line  relation  between  Alog  X  and  The  lower  portion 

to  the  table  lists  the  values  of  9,  the  excitation  temperatures  with  the 
probable  errors  derived  from  the  least  squares  calculations,  the  values 
of  log  u,  the  weighted  shifts,  the  abundances,  and  AlogT\. 


TABLE  Vi 
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Figure  9  illustrates  the  data  used  in  the  determination  of  the 
excitation  temperatures.  The  straight  lines  are  those  which  were  deter¬ 
mined  by  the  method  of  least  squares.  In  these  plots  the  sizes  of  the 
points  indicate  the  weights  they  were  given  in  the  calculations.  Circled 
points  indicate  zero  weights. 

The  curves  of  growth  obtained  for  the  four  models  are  shown  in 
Figure  10,  11,  12,  and  13.  The  theoretical  curves  are  represented  by 
the  solid  lines.  For  M->E  and  S-S  pure  scattering  the  theoretical  curves 
are  those  calculated  for  B°/B^  and  I°/I^  equal  to  2/3.  In  each  case, 
the  value  of  log  a  corresponding  to  the  pictured  theoretical  curve  is 
indicated  on  the  figure.  Also  indicated  is  the  value  of  log  c/v  used 
in  the  analysis,  In  .the  figures,  different  symbols  are  used  to  distin¬ 
guish  points  representing  observed  data  pertaining  to  different  atomic 
terms.  Following  the  procedures  outlined  in  this  study,  each  such  set 
of  points  was  fitted  individually  on  the  theoretical  curve. 


Figure  9.  Excitatlqn  Temperatures  Derived  prom  Ca  I  Lines 
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Results  for  T1  1 


Ti  I^NBS 

The  data  pertaining  to  the  Ti  I-NBS  lines  are  given  in  Table  VII. 
The  values  of  log  c/v  were  determined  from  the  a^F  plot,  which  fits 
around  the  knee  of  the  curve  of  growth.  No  lines  strong  enough  to 
deitermine  log  a  were  available,  and  the  values  listed  in  the  second 
row  of  Table  VII  merely  indicate  the  theoretical  curves  used  for  fit¬ 
ting  the  observed  plots.  To  call  attention  to  this,  the  values  of  log 
a  listed  in  the  table  are  enclosed  in  parentheses.  A  value  of  Alog  7^ 
for  titanium  is  not  given  since  the  abundances  obtained  for  Ti  II  also 
represent,  as  discussed  earlier,  the  total  abundances  for  titanium. 

Figure  14  illustrates  the  excitation  temperature  data,  and  Figures 
J.5,  16,  17,  and  18  are  the  curves  of  growth  obtained  for  Ti  I-NBS. 
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Figure  14.  Excitation  Temperatures  Derived  From  Ti  I-NBS  Lines. 
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Figure  16.  Milne-Eddington  Pure  Absorption  Curve  of  Growth  for  Ti  I-NBS, 
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Figure  18.  Schuster- Schwarzschild  Pure  Absorption  Curve  of  Growth  for  Ti  I-NBS. 
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T1  I- King 

The  Ti  I- King  data  are  presented  in  Table  VIII.  The  values  of 
log  c/v  were  determined  from  the  a^F  plot.  Values  of  log  a  could  not 
be  determined  because  of  the  absence  of  strong  lines,  and  the  values 
listed  indicate  the  theoretical  curves  used  for  fitting  the  observed 
plots . 

Since  the  f-values  of  the  Kings  are  on  a  relative  scale  it  was 
necessary  to  determine  a  value  of  Alog  gf  for  use  in  the  abundance 
calculations.  The  absolute  f-values  of  the  NBS  were  not  used  to  obtain 
a  value  of  Alog  gf  since  the  accur  -y  of  the  NBS  values  is  open  to  some 
question,  and  it  was  desired  to  keep  the  determinations  for  Ti  I-NBS 
and  Ti  I- King  as  independent  as  possible.  Allen's  (1955)  absolute 
f-values  were  used  for  this  purpose,  giving  a  value  of  Alog  gf  equal 
to  -3.30.  The  absolute  f-values  listed  by  Allen  were  determined  from 
the  Kings'  relative  values  using  f-sum  rules  and,  due  to  the  scarcity 
of  accurate  measurements  of  absolute  f-values  for  Ti  I,  were  regarded 
as  the  best  available. 

The  plots  of  the  data  used  to  determine  the  excitation  temperatures 
are  presented  in  Figure  19,  and  the  curves  of  growth  are  illustrated 
in  Figures  20,  21,  22,  and  23. 
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Figure  19.  Excitation  Temperatures  Derived  From  Ti  1-King  Lines. 
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Figure  23.  Schuster- Schwarzschild  Pure  Absorption  Curve  of  Growth  for  Ti  I- King. 
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Results  for  Ti  II 

NBS  f-values  and  values  derived  from  Wright's  solar  log  X^'s  were 

combined  for  use  in  the  study  of  the  Ti  II  lines.  Table  IX  gives  the 

_  2 

results  obtained.  Since  the  values  of  v  were  the  same  for  the  a  P  and 

^e 

4 

b  P  terms,  these  were  combined  to  form  a  single  plot  defining  a  large 

portion  of  the  cnrve  of  growth.  This  plot  was  used  to  determine  both 

2  2 

the  values  of  log  c/v  and  log  a.  The  b  D  and  a  H  terms  were  also  com¬ 
bined  since  the  respective  values  of  were  practically  the  same. 

The  excitation  temperatures  found  from  this  set  of  lines  are  very 
uncertain.  This  is  made  apparent  by  the  four  very  different  values 
obtained  from  the  application  of  the  four  idealized  models.  The  range 
of  represented  by  the  Ti  II  lines  is  quite  small,  In  such  a  situation, 
small  changes  in  the  values  of  Alog  X  can  result  in  a  relatively  large 
change  in  the  slope  of  the  best  straight  line  relating  Alog  X  and 
thus  producing  a  large  variation  in  the  excitation  temperature. 

Figure  24  gives  the  plots  of  the  data  pertaining  to  the  excitation 
temperature  determinations.  Figures  25,  26,  27,  and  28  are  the  curves 
of  growth  obtained  for  Ti  II. 
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Figure  24.  Excitation  Temperatures  Derived  From  Ti  II  Lines. 
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Figure  26.  Milne-Eddington  Pure  Absorption  Curve  of  Growth  for  Ti  II. 
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Figure  27.  Schuster- Schwarzschild  Pure  Scattering  Curve  of  Growth  for  Ti  II. 
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Figure  28.  Schuster-Schwarzschild  Pure  Absorption  Curve  of  Growth  for  Ti  II. 


97 


Results  for  Cr  I 

Cr  I-NBS 

The  results  for  Cr  I-NBS  are  given  In  Table  X.  The  a^D  plot  was 
used  to  determine  log  c/v.  The  scatter  of  the  points  In  the  a^D  plot 
was  quite  large,  making  the  determination  of  log  c/v  rather  uncertain. 
Also,  the  points  representing  the  two  weakest  lines  of  the  a^D  term 
were  separated  a  large  distance  from  the  other  points  on  the  plot. 

Hence,  these  two  weak  lines  exerted  considerable  Influence  In  the  deter¬ 
mination  of  log  c/v.  Since  the  possible  effects  of  blending  are  most 
serious  for  weak  lines,  this  influence  leads  to  uncertainty  in  the 
values  of  log  c/v. 

An  attempt  was  made  to  determine  the  values  of  log  a  by  comparing 
plots  of  Alog  X  versus  Xg  from  fittings  on  theoretical  curves  with 
various  values  of  log  a  to  see  which  gave  the  best  straight  lines.  The 
values  of  log  a  decided  upon  in  this  manner,  and  listed  in  Table  X,  are 
very  uncertain  due  to  the  large  scatter  of  the  points. 

Figures  29,  30,  31,  32,  and  33  present  the  plots  illustrating  the 
least  squares  determinations  of  the  excitation  temperatures  and  the 
curves  of  growth  obtained  for  each  of  the  four  models. 
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Figure  29.  Excitation  Temperatures  Derived  From  Cr  l-NBS  Lines. 
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Figure  30.  Milne-Eddington  Pure  Scattering  Curve  of  Growth  for  Cr  I-NBS, 


Figure  31.  Milne-Eddington  Pure  Absorption  Curve  of  Growth  for  Cr  I-NBS 
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Figure  32 . 
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Figure  33.  Schuster- Schwarzschild  Pure  Absorption  Curve  of  Growth  for  Cr  I-NBS 
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Cr 

The  Cr  I-Hill  data  are  presented  in  Table  XI.  The  values  of  log  c/v 
were  determined  from  the  a^D  plot.  The  values  of  log  a  were  determined 
in  the  same  manner  as  for  Cr  I-NBS,  and  are  also  quite  uncertain. 

Hill's  f- values  are  on  a  relative  scale,  and  to  obtain  a  value  of 
Alog  gf  for  use  in  the  abundance  calculations  the  absolute  f-values  of 
Bell  et  al.  (1958a)  were  employed.  These  measurements  are  regarded  as 
highly  accurate.  The  value  obtained  for  Alog  gf  was  -3.42. 

Figure  34  presents  the  plots  of  the  data  and  least  squares  straight 
lines  used  in  finding  the  excitation  temperatures.  The  curves  of  growth 
for  Cr  I-Hill  are  reproduced  in  Figures  35,  36,  37,  and  38. 


106 


X 

Cd 

iJ 

M 

< 

H 


w 


h4 

H 

X 

I 


u 

o 


X 

X 

Q 

Pd 

> 

H 

X 

Pd 

P 

I 


b 

O 


'0) 

ro 

:0 

w 


(U 

•a 

I 

X 


u 

4-1 

CNJ 

QO 

1-^ 

10 

CM 

00 

CM 

vO 

in 

SO 

1—4 

c 

•H 

00 

ON 

ON 

0 

P^ 

VO 

so 

0 

m 

o 

X 

« 

• 

• 

P 

• 

• 

• 

• 

• 

• 

• 

•H 

CO 

CO 

CO 

CO 

CM 

CO 

CO 

n 

n 

CO 

CO 

U 

m 

a 

o 

00 

M 

• 

• 

o 

lO 

rH 

X! 

(0 

1 

cs 

0 

n 

ON 

in 

00 

0 

0 

p^ 

m 

X) 

bO 

00 

ON 

in 

0 

0 

<!• 

CM 

CM 

r-l 

F— 1 

< 

O 

• 

• 

• 

• 

■ 

• 

• 

• 

• 

• 

n 

CM 

CM 

0 

0 

d 

0 

0 

a 

0 

0 

< 

1 

1 

u 

00 

CO 

0 

m 

CJO 

C30 

CM 

0 

^1 

0 

bO 

•H 

vO 

SO 

ON 

ON 

NO 

vO 

NO 

NO 

ON 

c 

JC 

• 

• 

• 

• 

• 

• 

• 

• 

• 

•H 

CO 

CO 

CO 

tn 

CM 

CM 

CO 

CO 

CO 

CO 

en 

CO 

o 

<1) 

00 

u 

• 

• 

4J 

m 

X 

CQ 

t 

00 

ON 

1—4 

in 

cn 

in 

^4 

so 

G 

bO 

vD 

nD 

'O 

1—4 

0 

CO 

CO 

CM 

fH 

'j- 

1-^ 

QO 

O 

• 

« 

• 

• 

• 

• 

• 

* 

• 

f 

* 

rH 

CO 

CM 

CM 

0 

0 

0 

0 

0 

0 

0 

0 

< 

1 

1 

U 

4-1 

00 

vt 

0 

CM 

in 

in 

00 

vO 

ON 

c 

•H 

0 

00 

ON 

0 

00 

p^ 

P>. 

vO 

f-4 

in 

0 

X 

• 

• 

• 

• 

• 

■ 

• 

• 

• 

• 

• 

CO 

CO 

<!■ 

tn 

CM 

CO 

CO 

CO 

cn 

CO 

<)■ 

CO 

o 

a 

o 

00 

t 

• 

o 

lO 

X 

CO 

1 

00 

in 

m 

CO 

CM 

ON 

00 

CM 

0 

<1* 

0 

43 

bO 

ON 

s£> 

0 

1—4 

>0 

CM 

CM 

m 

r-4 

<: 

O 

• 

• 

• 

• 

» 

• 

» 

• 

• 

• 

• 

1—4 

CO 

CM 

CM 

0 

0 

0 

0 

0 

0 

0 

0 

< 

1 

1 

u 

4-1 

v£> 

o^ 

in 

0 

0 

SO 

F— 1 

p^ 

SO 

bO 

•H 

r>- 

00 

0 

0 

00 

00 

P^ 

p^ 

0 

C 

X 

• 

• 

• 

• 

• 

• 

• 

« 

• 

• 

• 

•H 

CO 

CO 

CO 

cn 

CO 

CO 

CO 

CO 

en 

CO 

CO 

M 

in 

(U 

o 

00 

4J 

♦ 

• 

U 

in 

XJ 

CQ 

1 

so 

0 

rv 

00 

1“^ 

CO 

CM 

so 

CM 

m 

CO 

0 

bO 

r*- 

00 

v£> 

r— 4 

0 

in 

'i- 

CO 

CM 

m 

0 

C/3 

o 

♦ 

• 

• 

* 

• 

• 

• 

• 

• 

• 

• 

CO 

CM 

CM 

0 

0 

0 

0 

0 

0 

0 

0 

< 

1 

4: 

bO 

•H 

1—4 

fO 

0 

0 

CM 

0 

F^ 

F^ 

r-4 

G 

3 

44 

•  > 

0  CO 

<i» 

u 

i  ^ 

•  a 

fO 

CM 

in 

1—4 

CM 

CO 

CM 

1—1 

en 

CM 

CM 

0  -H 

^4 

bO 

bO 

Z  P 

o 

o 

0 

o^ 

tn 

0 

ON 

00 

Q 

00 

fH 

n- 

<u 

0 

ON 

o> 

in 

pv 

oo 

<Ti 

5 

0 

f-4 

IX 

• 

• 

• 

• 

• 

• 

• 

• 

0 

0 

0 

CM 

CM 

nj 

CM 

CO 

en 

CO 

cn 

e 

0 

0 

B 

CO 

(O 

Q 

0 

X 

X 

J® 

p 

0 

b 

Pd 

G 

in 

in 

m 

m 

ps. 

n 

in 

cn 

cn 

P^ 

H 

(0 

(0 

(D 

(0 

(0 

<0 

X 

CO 

CO 

N 

m 

m 


CM 


cn 


CO 

CO 


CD 


Excitation  Temperature  (°K)  4449  ±  146  4335  ±  182  4484  ±  151  4370  ±  182 


108 


Xg 

Figure  34.  Excitation  Temperatures  Derived  From  Cr  I-Hill  Lines. 


Milne-Eddington  Pure  Scattering  Curve  of  Growth  for  Cr  I-Hill 


Ill 


Figure  37.  Schuster- Schwarzschild  Pure  Scattering  Curve  of  Growth  for  Cr  I-Hill, 
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Results  for  Mn  1 

The  NBS  f-values  were  employed  for  Mn  I.  Table  XII  gives  the 
results  obtained  from  the  curve  of  growth  analysis .  The  values  of 
log  c/v  were  adopted  from  those  determined  for  Fe  I-NBS.  The  values 
of  log  a  were  taken  to  be  those  that  gave  the  best  straight  line  rela¬ 
tions  between  Alog  X  and 

The  plots  illustrating  the  determinations  of  the  excitation 
temperatures  are  shown  in  Figure  39.  The  curves  of  growth  are  shown 
in  Figures  40,  41,  42,  and  43. 
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Figure  39.  Excitation  Temperatures  Derived  From  Mn  I  Lines. 
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Figure  40.  Milne- Eddington  Pure  Scattering  Curve  of  Growth  for 
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Pure  Absorption  Curve  of  Growth  for 
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Results  for  Fe  I 


Fe  I^NBS 

The  results  for  Fe  I-NBS  are  presented  in  Table  XIII.  The  values 

5  3 

of  log  c/v  and  log  a  were  determined  from  the  a  P  and  a  F  plots.  The 

5  3 

a  P  plot  fit  around  the  knee  of  the  curve  of  growth  and  the  a  F  plot 

determined  the  damping  and  transition  portions  of  the  curve. 

Figures  44,  45,  46,  47,  and  48  give  the  excitation  temperature 

plots  and  the  curves  of  growth. 
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Figure  45.  Milne-Eddington  Pure  Scattering  Curve  of  Growth  for  Fe  I-NBS 
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Figure  47.  Schuster-Schwarzschild  Pure  Scattering  Curve  of  Growth  for  Fe  I-NBS. 


Figure  48.  Schuster- Schwarschild  Pure  Absorption  Curve  of  Growth  for  Fe  I-NBS 
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I-Klne 

Table  XIV  contains  the  results  obtained  for  Fe  I- King.  The  values 

5  3 

of  log  c/v  and  log  a  were  determined  from  the  a  D  and  a  F  plots. 

The  reliable  absolute  f-values  of  Bell  et  al .  (1958b)  were  used  to 
obtain  a  value  of  Alog  gf  equal  to  >3.27  for  use  In  the  abundance 
determinations , 

Excitation  temperature  data  are  plotted  In  Figure  49.  The  curves 
of  growth  are  shown  In  Figures  5D,  51,  52,  and  53. 
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Growth  for  Fe  1-King. 
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Figure  51.  Milne-Eddington  Pure  Absorption  Curve  of  Grovth  for  Fe  I- King. 
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ure Absorption  Curve  of  Growth  for  Fe  I- King. 
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Fe  I-Carter 

Table  XV  presents  the  results  obtained  for  Fe  I-Carter.  The  values 
of  log  c/v  were  taken  to  be  the  same  as  the  values  obtained  for  Fe  I- 
King.  Values  of  log  a  could  not  be  determined  because  of  the  absence 
of  strong  lines.  In  Table  XV,  the  listed  values  of  log  a  indicate  the 
theoretical  curves  used  in  the  analysis. 

Carter's  f- values  were  intended  to  be  on  an  absolute  scale.  How¬ 
ever,  the  absolute  scale  used  by  Carter  is  now  believed  to  be  incorrect, 
and  for  the  purposes  of  this  investigation  the  scale  was  regarded  as 
relative.  Carter  first  determined  his  f-values  on  a  relative  scale, 
that  of  the  Kings  (1938),  then  converted  to  the  absolute  scale  of  R,  B. 
King  (1942).  Later,  as  discussed  in  the  paper  by  Bell  et  al.  (1958b), 
it  was  necessary  to  revise  King's  absolute  scale  because  of  improved 
vapor  pressure  data.  Due  to  this  revision,  a  correction  is  required 
in  the  absolute  scale  used  by  Carter.  This  was  taken  into  account 
in  the  calculation  of  abundances  by  applying  a  value  of  Alog  gf  equal 
to  0.46. 

The  excitation  temperature  plots  and  the  curves  of  growth  are 
shown  in  Figures  54,  55,  56,  57,  and  58. 
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Figure  54.  Excitation  Teinperaturea  Derived  From  Fe  I-Carter  Lines. 
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tion  Curve  of  Growth  for  Fe  1-Carter. 
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Scattering  Curve  of  Growth  for  Fe  I-Carter. 
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Results  for  Co  I 

For  Co  J  the  NBS  f-values  were  us'^d.  Table  XVJ  presents  the  results 
obteined,  The  valutes  of  log  c/v  were  adopted  from  those  obtained  for 
Fe  X-INBS,  Because  of  the  absence  of  strong  lines  it  was  not  possible 
to  determine  values  of  log  a,  The  values  listed  in  Table  XVI  merely 
indicate  the  theoretical  curves  used  in  the  analysis. 

The  excitation  temperature  data  are  plotted  in  Figure  59,  The 
curves  of  growth  are  shown  in  Figures  60,  61,  62,  and  63, 


TABLE  XVI 
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Figure;  59.  Excitation  Temperatures  Derived  From  Co  I  Lines. 
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log  Tl^ 

Figure  60.  Milne- Eddington  Pure  Scattering  Curve  of  Growth  for  Co 


Figure  63.  Schuster- Schvarzschild  Pure  Absorption  Curve  of  Growth  for  Co 
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Results  for  Ri  I 

The  fwvalues  used  here  were  those  of  the  NBS.  Table  XVII  gives 
the  results  obtained  for  Ni  I,  The  vslues  of  log  c/v  listed  in  the 
table  vere  adopted  from  the  values  obtained  for  Fe  I-'NBS.  Because  of 
the  pbsence  of  strong  lines,  values  of  log  a  could  not  be  determined, 
end  the  values  listed  in  Table  XVII  merely  refer  to  the  theoretical 
curves  used  in  the  analysis* 

figures  6A,.  65,  66,  67,  and  68  present  tbe  excitation  temperature 
plots  and  the  curves  of  growth  obtained  for  each  of  the  four  models. 
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Xe 


Flgujre  64.  Excitation  Temperatures  Derived  From  Ni  I  Lines. 
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gure  65.  Milne- Eddington  Pure  Scattering  Curve  of  Growth  for  Ni 


Figure 
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Figure  67.  Schuster- Schwarzschild  Pure  Scattering  Curve  of  Growth  for  Ni 
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Figure  68.  Schuster- Schwarzschild  Pure  Absorption  Curve  of  Growth  for  Ni 


CHAPTER  V 


SUMMARY  AND  CONCLUSIONS 

The  results  of  this  investigation  are  summarized  in  Table  XVIII. 
Values  of  log  c/v,  or  log  a,  or  both  could  not  be  determined  from  the 
observed  curves  of  growth  for  some  of  the  sets  of  lines.  In  such  cases, 
values  of  log  c/v  and  log  a  were  adopted,  as  discussed  in  Chapter  IV, 
for  use  in  the  analysis.  In  Table  XVIII  such  values  are  enclosed  in 
parentheses  to  indicate  that  they  are  not  definitive  results  of  this 
investigation.  For  titanium  no  value  of  Alog  Tl  is  given  since  the 
abundances  listed  for  Ti  II  also  represent,  to  the  limit  of  the  accuracy 
attained  in  this  study,  the  total  abundances  obtained  for  titanium. 

The  results  obtained  from  the  Fe  I-NBS  lines  are  regarded  as  the 
most  reliable;  more  lines  were  available  here  than  for  any  other  set, 
and  the  observed  curves  of  growth  for  Fe  I-NBS  were  very  well  defined. 

A  curve  of  growth  analysis  is  a  statistical  procedure,  and  the  Fe  I-NBS 
data  must  be  considered  as  best  from  a  statistical  point  of  view. 

The  excitation  temperatures  found  from  the  Ti  II  lines  are  very 
uncertain,  as  evidenced  by  the  four  very  different  values  obtained  from 
the  application  of  the  four  idealized  models.  As  discussed  in  Chapter  IV 
in  the  presentation  of  the  results  for  Ti  II,  the  responsibility  for 
bhis  uncertainty  lies  in  the  small  range  of  represented  by  the  Ti  II 
lines.  In  such  a  situation,  small  changes  in  the  values  of  Alog  X  can 
result  in  a  relatively  large  change  in  the  slope  of  the  best  straight 
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TABLE  XVIII 
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line  relating  Alog  X  and  Xgi  thus  producing  a  large  variation  in  the 
excitation  temperature. 

Comparison  of  Results  Obtained  Using  Different  Sets  of  f-Values 

Most  of  the  lines  included  in  the  Ti  I- King,  Cr  I- Hill,  Fe  I-King, 

3nd  Fe  I-Carter  data  were  also  present  among  the  NBS  lines.  For  these 
common  lines  the  NBS  f-values  are  compared  in  Figure  69  with  those  from 
the  other  sources.  The  NBS  f-values  are  on  an  absolute  scale,  while 
those  from  the  other  sources  are  on  relative  scales.  The  values  of 
Alog  gf  used  in  the  abundance  calculations  for  Ti  I-King,  Cr  I-Hill, 

Fe  I-King,  gnd  Fe  I-^Carter  are  the  quantities  necessary  to  convert  from 
the  relative  to  absolute  scales.  This  can  be  seen  from  Equation  (31) 
of  Chapter  IV.  However,  there  is  no  guarantee  that  the  absolute  scales 
determined  in  this  manner  will  be  the  same  as  the  NBS  scale.  In  Figure  69 
the  values  of  Alog  gf  were  used  to  establish  the  positions  of  the  45- 
degree  lines.  Points  falling  exactly  on  the  lines  represent  perfect 
agreement  between  the  absolute  f-values  of  the  NBS  and  the  derived 
absolute  f-values  of  the  other  sources.  Some  large  deviations  from  agree¬ 
ment  are  apparent  in  the  figure.  At  the  appropriate  points  in  the  follow¬ 
ing  discussion  these  deviations  are  taken  into  consideration. 

For  Ti  I,  Cr  I,  and  Fe  I  many  more  lines  were  available  with  NBS 
f-^values  than  with  f-values  from  the  other  sources;  hence  small  differ¬ 
ences  might  be  expected  in  the  results  obtained  using  the  different  sets 
of  f-values.  Falling  into  this  category  are  variations  in  the  values  of 
log  c/v  and  the  excitation  temperatures  found  from  the  Ti  I-NBS  and  Ti  I- 
King  lines  and  in  the  values  of  log  c/v  found  from  the  Fe  I-NBS  and  Fe  I- 
King  lines.  However,  certain  differences  are  present  which  cannot  be 


log  gfX  (NBS)  log  gfX  (NBS) 
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Figure  69.  Comparison  of  the  NBS  f- Values  With  (a)  the  Ti  I- King  f- Values 
(b)  the  Cr  I-Hill  f-Values,  (c)  the  Fe  I-King  f-Values,  and  (d)  the  Fe  I 
Carter  f-'Values  for  the  Lines  Used  in  This  Investigation. 
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attributed  simply  to  statistical  causes.  These  are  discussed  below. 

Ti  I 

The  most  important  differences  between  the  Ti  I-NBS  and  the  Ti  I- 
Kipg  results  are  in  the  abundances.  As  is  evident  in  Figure  69a,  the 
absolute  scale  of  the  NBS  f-values  does  not  agree  with  the  absolute  scale 
derived  for  the  Ti  I-King  f-values  from  Alog  gf  =  -3.30.  This  disagree¬ 
ment  may  be  largely  responsible  for  the  apparent  differences  in  the  abun^ 
dances.  The  fault  may  lie  in  the  value  of  Alog  gf,  in  the  absolute  scale 
used  by  the  NBS,  or  in  a  combination  of  the  two. 

Cr  I 


The  Cr  I-NBS  and  Cr  I-Hill  results  differ  markedly  in  the  values  of 
log  c/v,  the  excitation  temperatures,  and  the  abundances. 

As  remarked  in  Chapter  IV  in  the  presentation  of  the  results  obtained 
from  the  Cr  I-NBS  lines,  two  weak  lines  were  very  influential  in  the 
determinations  of  log  c/v.  These  two  lines  were  not  included  among  the 
Cr  I-Hill  lines.  The  large  differences  between  the  Cr  I-NBS  and  the 
Cr  I-Hill  values  of  log  c/v  aye  due  to  the  effects  of  these  two  lines. 

The  differences  between  the  Cr  I-NBS  and  Cr  I-^Hill  excitation 
temperatures  are  a  result  of  the  non-uniform  relationship  between  the  f- 
values  used  in  the  two  cases.  It  is  apparent  from  Figure  69b  that  the 
f-- values  for  Cr  I  fall  into  two  groups:  those  for  lines  with  Xg  greater 
than  2.53  ev  and  those  for  lines  with  Xg  less  than  1.03  ev.  The  excitation 
temperature  is  found  from  the  slope  of  the  best  straight  line  relating 
Alog  X  and  x^.  The  values  of  Alog  X  depend  upon  the  f-values,  e.g., 

Alog  X  “  log  -  log  gf\  in  the  S-S  pure  scattering  model.  For 
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convenience  in  the  discussion,  consider  thet  the  Cr  I-Hill  f-values  are 
placed  on  the  absolute  scale  given  by  Alog  gf  =  -3.42.  (Whether  the 
scale  is  relative  or  absolute  does  not  affect  the  excitation  temperature 
determinations,  but  the  uniformity  of  the  scale  is  important.)  Then,  for 
lines  with  less  than  1.03  ev,  the  values  of  Alog  X  obtained  using  the 
Cr  I-NBS  f-values  should  be  approximately  the  same  as  those  obtained 
using  the  Cr  I-Hill  f-values.*  This  can  be  seen  from  Figure  69b.  For 
Xg  greater  than  2.53  ev,  the  values  of  Alog  X  obtained  using  the  Cr  IrNBS 
f-values  should  be  considerably  larger  than  those  obtained  from  the 
Cr  I-Hill  f-values.  Hence,  the  slope  of  the  best  straight  line  relating 
Alog  X  and  x^  obtained  using  the  Cr  I-NBS  f-values  should  be  substantially 
different  from  that  obtained  using  the  Cr  I-Hill  f-values.  This  is  a 
possible  explanation  of  the  differences  between  the  excitation  temperatures 
found  from  the  Cr  I-NBS  and  the  Cr  I-Hill  lines. 

The  abundances  determined  from  the  Cr  I-NBS  lines  are  smaller  than 
those  found  from  the  Cr  I-Hill  lines.  The  abundances  in  each  of  the  four 
models  are  dependent  upon  the  quantity 

shift  +  log  u  -  log  c/v. 

The  shift  is  given  by  Alog  X  +  (5040/T)Xg,  where  T  is  the  excitation  tem¬ 
perature.  Larger  excitation  temperatures  give  smaller  shifts,  and  thus 
smaller  abundances.  This  effect  is  partially  offset  by  the  increase  of 
log  u  with  temperature.  However,  the  variation  of  log  u  with  temperature 

*This  assumes  that  the  same  vertical  shift  (log  c/v)  is  used  for  both 
cases.  Although  this  was  not  true  in  the  actual  analysis  of  the  Cr  I 
lines,  inclusion  of  this  consideration  would  only  complicate  the  argument 
and  not  alter  the  outcome.  This  can  be  seen  by  noting  that  a  change  in 
log  c/v  is  accompanied  by  a  change  in  Alog  X,  and  for  Cr  I,  where  most  of 
the  lines  are  weak,  most  of  the  changes  in  Alog  X  should  be  approximately 
the  same.  Hence,  the  slope  of  the  best  straight  line  relating  Alog  X 

and  X  should  not  be  altered  appreciably. 

e 
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is  comparatively  small.  Hence,  the  much  greater  excitation  temperatures 
obtained  for  Cr  I-NBS  than  for  Cr  I-Hill  produced  to  a  large  extent  the 
differences  in  the  calculated  aburjdances. 

Another  consideration  of  importance  in  the  abundance  calculations  is 
the  connection  between  Alog  X  and  log  c/v.  By  considering  the  shape  of 
the  curve  of  growth,  it  can  be  seen  that  as  log  c/v  increases  so  does 
Alog  X.  In  the  region  well  below  the  knee  of  the  curve  of  growth,  the 
slope  of  the  curve  is  unity.  For  lines  confined  to  this  region,  increases 
in  log  c/v  are  accompanied  by  like  increases  in  Alog  X,  and  the  effects 
of  these  changes  exactly  cancel  one  another  in  the  abundance  calculations. 
For  stronger  lines,  hofwever,  an  increase  in  log  c/v  produces  a  larger 
increase  in  Alog  X,  with  a  resulting  increase  in  the  calculated  abundance. 
In  the  case  of  Cr  I  most  of  the  lines  were  weak,  and  it  would  be  expected 
that  the  differences  between  the  values  of  log  c/v  found  for  Cr  I~NBS 
and  Cr  I-Hill  would  be  largely  compensated  in  the  abundance  calculation's 
by  the  resulting  differences  in  the  values  of  Alog  X.  This  would  be  true 
were  it  not  for  the  non- uniformity  in  the  relationship  between  the  f-values 
used  for  the  two  sets  of  lines,  and  the  effect  of  this  on  the  values  of 
Alog  X.  Because  of  this  non- uniformity  the  values  of  Alog  X  found  from 
the  Cr  I-NBS  lines  with  Xg  greater  than  2.53  ev  were  disproportionately 
larger  than  the  corresponding  Cr  I-Hill  values.  This  acted  in  such  a 
way  as  to  decrease,  but  not  nullify,  the  effect  in  the  abundance  calcu¬ 
lations  of  the  differences  between  the  excitation  temperatures. 

Fe  1 


Important  differences  exist  in  the  excitation  temperatures  found 
from  the  three  sets  of  Fe  I  lines.  These  differences  seem  to  reflect 
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the  differences  in  the  ranges  of  the  excitation  potentials  applying  in 
each  case.  For  Fe  I-King  the  values  of  Xg  ranged  from  0.05  to  1.54  ev, 
for  Fe  I-Carter  from  0.11  to  3.56  ev,  and  for  Fe  I-NBS  from  0.03  to  4.53 
ev.  As  a  point  of  interest,  excitation  temperatures  were  calculated  from 
the  Fe  I-NBS  values  of  Alog  X,  obtained  from  the  S-S  pure  scattering 
model,  covering  first  the  same  range  in  Xg  as  represented  by  the  Fe  I-' 

King  lines,  and  second  the  same  range  in  Xg  as  represented  by  the  Fe  I- 
Carter  lines.  For  Xg  ia  the  range  0.05  to  1.54  ev  the  temperature  found 
was  3616°K,  and  for  Xg  the  range  0.11  to  3.56  ev  the  temperature  deter¬ 
mined  was  5009°K.  The  first  value  corresponds  closely  to  that  obtained 
for  the  Fe  I-King  temperature  (3732°’''  for  S-S  pure  scattering),  and  the 
second  is  not  greatly  different  from  the  Fe  I-Carter  value  (4578°K) . 

It  appears  that  the  value  determined  for  the  excitation  temperature 
depends  upon  the  range  of  Xg>  with  lower  temperatures  being  characteristic 
of  the  lower  levels  of  excitation.  This  lends  weight  to  the  reasonable 
assumption  that  lines  arising  from  the  higher  levels  in  the  atom  are 
more  likely  to  be  formed  principally  in  deeper  layers  of  the  atmosphere 
where  the  temperature  is  higher.  In  a  curve  of  growth  analysis  the 
assumption  is  made  that  a  single  value  of  the  excitation  temperature  can 
be  specified  which  yields  the  correct  distribution  of  atoms  among  the 
various  energy  levels.  As  is  apparent  from  the  above  discussion,  this 
assumption  may  be  open  to  serious  question  when  a  wide  range  of  excitation 
potentials  is  involved. 

As  described  in  the  discussion  of  the  Cr  I  results,  differences  in 
the  excitation  temperatures  have  a  rather  strong  influence  in  producing 
differences  in  the  calculated  abundances.  However,  this  is  not  apparent 
in  the  abundances  found  from  the  three  sets  of  Fe  I  lines.  In  the 
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abundance  calculations  for  Fe  I-N3S  and  Fe  I-King,  the  differences  in  the 
values  of  log  c/v  exert  a  greater  influence  than  the  differences  in 
excitation  temperature.  Most  of  the  lines  included  in  these  two  sets 
are  strong  lines,  and  for  strong  lines  an  increase  in  log  c/v  produces 
a  larger  increase  in  Alog  X.  Since  the  Fe  I-NBS  values  of  log  c/v  are 
larger  than  the  Fe  1-King  values,  this  effect. acts  in  a  manner  contrary 
to  that  of  the  excitation  temperatures  in  the  abundance  calculations  for 
Fe  I-NBS  and  Fe  I-King,  Also,  the  fact  that  smaller  values  of  log  a 
were  obtained  from  the  Fe  I-NBS  lines  than  from  the  Fe  I-King  lines* 
(except  for  the  S-S  pure  scattering  case)  tends  to  cause  further  increases 
in  the  Fe  I-NBS  values  of  Alog  X  over  those  of  Fe  I-King.  The  result  is 
that  the  Fe  I-NBS  abundances  are  larger,  except  for  the  S-S  pure  scat¬ 
tering  case,  than  the  Fe  I-King  abundances. 

In  Figure  69d  it  is. apparent  that  the  derived  absolute  scale  for 
the  Fe  I-Carter  f-values  does  not  agree  with  the  absolute  scale  of  the 
Fe  1-NBS  f-values.  This  disagreement  is  mainly  responsible  for  the 
differences  between  the.  Fe  I-Carter  and  Fe  I-NBS  abundances . 

Comparison  of  Results  Obtained  Using  Different  Idealised  Models 

A  comparison  of  the  results  found  from  the  four  different  sets  of 
theoretical  curves  of  growth  reveals  well  defined  trends  in  the  values 
of  log  c/v  and  the  abundances. 

Except  for  Ti  II,  the  S-S  pure  scattering  model  produced  the  largest 
values  of  log  c/v,  and  the  M-E  pure  absorption  model  the  smallest  values. 

*The  different  values  of  log  a  were  brought  about  by  the  different 
values  of  log  c/v  obtained  from  the  two  sets  of  lines.  This  can  be  seen 
by  referring,  in  Chapter  IV,  to  the  description  of  the  procedures  used  in 
determining  log  a. 
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The  reason  for  this  is  apparent  from  an  examination  of  Figure  70,  where 
theoretical  curves  with  log  a  *  -1.8  for  the  four  models  are  compared. 

Pvpn  if  the  S-S  pure  scattering  curye  were  to  be  shifted  horizontally  so 
that  the  Doppler  portions  (i.e,,  the  portions  below  the  knees)  of  all  four 
curvep  were  coincident,  the  transition  portion  of  the  S-'S  pure  scattering 
curve  would  still  be  higher  than  the  transition  portions  of  the  other 
curves.  Hence,  the  largest  vertical  shifts,  and  thus  the  largest  values 
of  log  c/v,  should  result  from  the  application  of  the  S-S  pure  scattering 
curve.  Conversely,  the  smallest  value  should  be  given  by  the  M-E  pure 
absorption  curve. 

It  should  be  noted  that  for  use  with  the  M-E  pure  scattering  model, 

the  M-E  pure  absorption  model,  and  the  S~S  pure  absorption  model  the 

observed  curves  of  growth  were  altered  from  those  used  with  the  S-S  pure 

scattering  model  by  the  application  of  the  limb  darkening  corrections-- 

log  2R  and  Alog  Tl  .  However,  the  limb  darkening  corrections  in  most 
c  o 

cases  were  not  large  enough  to  change  the  trend  discussed  in  the  preceding 
paragraph. 

In  the  case  of  Ti  II  the  damping  parameter  was  important  in  the 
determination  of  log  c/v.  The  plot  used  to  determine  log  c/v  extended 
from  the  knee  well  into  the  damping  portion  of  the  curve  of  growth; 
therefore  this  plot  was  also  used  to  determine  log  a,  Log  a  =  -1.8  was 
found  for  both  scattering  models  and  log  a  =  -1.3  for  both  absorption 
models.  The  curves  for  log  a  =  -1.3  rise  more  rapidly  than  those  for 
log  a  ^  t1.8.  This,  combined  with  the  effects  produced  by  the  inclusion 
of  the  limb  darkening  corrections,  caused  the  S-S  pure  absorption  model 
to  give  the  largest  value  of  log  c/v  and  the  M-E  pure  scattering  model 


the  smallest  value. 
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The  values  of  log  N/p7t  found  from  the  M-E  models  are  larger,  in 
most  cases,  than  the  values  of  log  NH  determined  from  the  S-S  models. 

From  Figure  70  it  can  be  seen  that  the  curves  of  grovith  for  the  M-E 

models  tend  to  give  either  smaller  values  of  log  c/v  or  larger  values  of 

Alog  X  than  the  curve  for  the  S-S  pure  scattering  model.  The  result  in 
either  case  is  that  the  M-E  abundances  tend  to  be  larger  than  the  S-S 
pure  scattering  abundance.  The  S-S  pure  absorption  curves  of  Hunger  are 
plotted  with  log  3t^/2  as  abscissa  rather  than  simply  log  T^.  This 
amounts  to  shifting  the  curve  horizontally  so  that  the  Doppler  portion 
coincides  with  the  similar  portions  of  the  M-E  curves.  Hence,  the 
differences  in  the  values  of  Alog  X  obtained  from  this  model  and  the  M-E 
models  should  be  insignificant.  However,  the  term  log  3/2  appears  nega¬ 
tively  in  the  equation  giving  the  abundance  in  the  S-S  pure  absorption 
mpdel.  This  compensates  for  the  shift  in  the  abscissa,  and  the  resulting 
abundance  should  be  smaller  than  those  found  from  the  M-E  models. 

Compar:*  son  of  the  Results  of  This  Study  With  Those  of  Greenstein 

Greenstein  (1948),  in  his  study  of  a  group  of  F  stars,  employed 
observed  curves  of  growth  based  mainly  upon  Fe  I  lines  and  Wright's  solar 
log  Xj-values.  These  were  analyzed  using  theoretical  curves  of  growth 
based  upon  the  M-E  model.  Greenstein's  values  of  log  v,  and 

which  are  discussed  below,  were  obtained  from  his  observed  curve  of  growth 
for  Fe  I.  . 

Greenstein  found  for  9  Ursae  Majoris  a  value  of  log  v  =  5.44.  (The 
velocity  parameter  in  Greenstein's  notation  is  represented  by  "V.") 

This  leads  to  a  value  of  log  c/v  equal  to  5,04,  which  agrees  closely 
with  the  most  reliable  values  (the  Fe  1-NBS  values)  obtained  in  this 


investigation. 
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He  obtained  a  value  of  equal  to  1.7,  which  corresponds  to 

log  a  *  -2.3.  Here  F  is  the  effective  damping  constant  and  y^j^  is  the 
classical  damping  constant.  By  referring  to  Figure  8  of  Chapter  III  it 
can  be  seen  that  for  the  strongest  lines,  which  are  important  in  the 
determination  of  log  a,  the  values  of  log  V/X  measured  by  Greenstein 
are  smaller  than  those  measured  in  this  investigation,  This  difference 
could  have  caused  the  value  of  log  a  determined  by  Greenstein  to  be 
smaller  than  the  values  found  in  this  study, 

For  0  Greenstein  obtained  the  value  0.98,  which  corresponds  to 
exc 

an  excitation  temperature  equal  to  5150°K,  The  lower  level  excitation 
potentials  of  the  Fe  I  lines  used  by  Greenstein  range  from  0,00  to  4.09 
ev.  The  extent  of  this  range  is  intermediate  to  those  represented  by 
the  Fe  I-Carter  and  Fe  I-NBS  lines  used  in  this  investigation.  As 
expected  from  the  apparent  relationship,  discussed  earlier,  between  the 
excitation  temperature  and  the  range  of  Xg>  the  excitation  temperature 
obtained  by  Greenstein  falls  between  the  values  found  in  this  study  from 
the  Fe  I-Carter  and  Fe  I-NBS  lines. 

Greenstein  determined  abundances  in  the  stars  relative  to  those  in 
the  sun,  and  these  quantities  cannot  readily  be  compared  with  the  abun¬ 
dances  found  in  this  investigation.  The  abundances  obtained  by  Wright 
(1948)  for  the  star  Procyon,  which  has  a  spectrum  similar  to  that  of 
0  Ursae  Majoris,  are  of  the  same  order  of  magnitude  as  those  found  for 
6  Ursae  Majoris  in  this  study. 

Comparison  of  Wrubel's  and  Hunger's  M-E  Pure 
Absorption  Curves  of  Growth 


Hunger's  M-E  pure  absorption  curves  of  growth  were  employed  in  this 


173 


investigation  although  curves  for  the  same  model  by  Wrubel  were  available. 
Wrubel's  curves  were  calculated  using  Chandrasekhar's  exact  solution  of 
the  equation  of  transfer  and  therefore  are  regarded  as  the  best  available 
for  the  Mn’E  pure  absorption  model.  However,  calculations  were  made  for 
only  two  values  of  log  a,  -I  and  •?-3,  thus  limiting  the  usefulness  of 
these  curves. 

The  abscissa  of  both  Wrubel's  and  Hunger's  curves  is  log  T]^.  (In 

Hunger's  notation  the  abscissa  is  log  C,  but  for  M-E  pure  absorption  C  = 

T]^.)  The  ordinate  of  Wrubel's  curves  is  log  W/b,  while  for  Hunger's 

curves  the  ordinate  is  log  W/2R  b.  Wrubel  calculated  curves  for  two 

c 

values  of  2/3  and  10/3.  If  B°/B^  is  equal  to  2/3,  then  2R^  is 

equal  to  unity.  In  this  case  Hunger's  ordinate  is  simply  log  W/b,  the 
same  as  Wrubel's  ordinate.  By  interpolating  in  Hunger's  tables,  curves 
were  obtained  for  log  a  =  -1  and  log  a  •*  -3.  These  were  then  compared 
with  Wrubel's  curves  for  B°/B^  =  2/3.  The  calculated  points  of  Wrubel's 
and  Hunger's  curves  could  not  be  compared  individually  because  Wrubel 
and  Hunger  did  not  use  the  same  sets  of  values  of  T)^  in  their  calculations. 
To  compare  the  curves  for  the  same  value  of  log  a,  the  two  curves  were 
plotted,  using  the  same  scale,  and  one  superposed  upon  the  other.  There 
were  no  noticeable  differences  in  the  two  curves  with  log  a  =  -1  nor  in 
the  two  curves  with  log  a  =*  -3. 

The  ordinate  of  Hunger's  curves  can  be  written  as 

2^  -  log  ^  -  log  2R^.  (36) 

c 

If  B°/B^  is  equal  to  10/3,  then  log  2R^  is  equal  to  -0,477.  Hunger's 
ordinate  in  this  case  is 

log  -  log  ^  t  0.477.  (37) 

c 
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To  obtain  a  curve  which  relates  log  Tl^  and  log  W/b  for  B°/B^  =  10/3, 
it  i?  only  necessary  to  subtract  the  quantity  0,474- from  Hunger's 
tabulated  values  of  log  W/2R^b,  For  log  a  =>  -1  and  Ipg  a  =  -3  there  were 
np  noticeable  differences  in  the  curves  obtained  in  this  manner  and  the 
curves  calculated  by  Wrubel  for  B°/B^  =  10/3. 

The  above  discussion  does  not  eliminate  the  possibility  that  signi¬ 
ficant  differences  might  exist  between  the  M-E  pure  absorption  curves  of 
Hunger  end  those  calculated  using  Wrubel 's  method  for  values  of  log  a 
other  than  -^l  and  -3  and  values  pf  B°/B^  other  than  2/3  and  10/3,  However, 
the  suggestion  is  strong  that  for  any  values  of  log  a  and  B°/B^  the  curves 
should  be  substantially  the  same. 

The  four  sets  of  curves  of  growth  used  in  this  investigation  are 
based  upon  idealized  models  which  represent  extreme  approximations  as 
to  the  structure  of  the  stellar  atmosphere  and  the  mechanism  of  line 
formation.  None  of  these  models  can  be  expected  to  perfectly  characterize 
the  actual  structure  and  radiation  processes  of  the  atmosphere  of  @  Ursae 
Majoris . 

In  the  case  of  the  star  Procyon,  Schrpeder  (1958)  has  pointed  out 
that  the  S-S  model  is  to  be  preferred  for  the  lines  of  low  and  interme¬ 
diate  excitation  potential  since  these  lines  are  forrned  veil  above  the 
photosphere.  For  lines  of  higher  excitation  potential,  which  are  formed 
prlpiarlly  in  the  same  region  as  the  continuum,  Schroeder  concluded  that 
the  M*’E  model  would  be  more  appropriate.  Due  to  the  similarity  of 
Procyon  and  @  Ursae  Majoris,  it  seems  likely  that  the  above  remarks 
should  also  apply  to  9  Ursae  Majoris. 
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In  the  spectrum  of  9  Ursae  Majoris  the  central  intensities  of  the 
strongest  lines  are  not  zero.  This  indicates  that  the  lines  are  not 
formed  solely  by  the  scattering  process,  which  requires  zero  central 
intensities  for  the  strongest  lines.  The  lines  are  most  likely  formed 
by  a  combination  of  the  processes  of  scattering  and  absorption. 

Although  the  four  idealized  models  (i.e.,  M-E  and  S-S  pure  scattering, 
and  M-E  and  S-S  pure  absorption)  represent  extreme  approximations,  the 
pesults  obtained  in  this  investigation  from  the  application  of  the  four 
different  models  are  not  greatly  variant, 

In  abundance  determinations  the  quantity  most  desired  is  JJ,  the 
number  of  ato^s  per  cubic  centimeter,  N  is  inherent  in  the  expressions 
NH  and  N/pn,  which  are  the  "abundances"  determined  in  this  study.  H 
is  the  depth  of  the  reversing  layer  characterizing  the  S-S  model.  In 
the  M-E  model  l/pn  is  a  measure  of  the  geometrical  depth  of  the  atmos- 
phete  corresponding  to  infinite  optical  depth  in  the  continuum.  It  is 
clear  that  a  knowledge  of  H  or  l/pn  is  prerequisite  to  obtaining  a 
value  for  N.  H  depends  upon  excitation,  ionization,  and  varies  from 
element  to  element,  while  l/pn  may  vary  strongly  with  optical  depth,  At 
best,  average  values  of  H  and  l/pn  can  be  obtained.  Until  more  refined 
techniques  are  devised,  curve  of  growth  analyses  will  yield  only  NH  or 
N/pM>  bhe  number  of  atoms  in  a  one  square  centimeter  column  in  the  line 
of  sight • 

A  curve  of  growth  analysis  yields  information  of  a  statistical  nature 
concerning  the  properties  of  a  stellar  atmosphere,  Suph  information  pro¬ 
vides  a  basis  for  more  detailed  studies  of  the  individual  line  profiles 
using  model  atmosphere  techniques.  Studies  of  this  type  are  called  for 
in  the  further  investigation  of  the  atmosphere  of  0  Ursae  Majoris, 


REFERENCES 


Allen,  C.  W.  1955,  Astrophvslcal  Quantities  (London:  Athlone  Press), 

Aller,  L.  H.  1960,  Stellar  Atmospheres,  ed .  J,  L.  Greenstein 
(Chicago:  The  University  of  Chicago  Press),  pp.  199,  246. 

Bell,  G,  D-,  Davis,  M,  H.,  Kirig,  R,  B.,  and  Routly,  P,  M.  1958a, 

Trans ,  I .  A .  U . .  10.  221. 

_ _ ,  1958b,  Ap.  J..  127.  775. 

Carter,  W,  W.  1949,  Phvs .  Rev . .  76.  962. 

Corliss,  C.  H.,  and  Eozmajt,  W.  R.  1962,  Experimental  Transition  Prob¬ 
abilities  for  Spectral  Lines  of  Seventy  Elements  (NBS  Mono,  No.  53 
[Washington:  Government  Printing  Office]), 

Deutsch,  A.  J.  1954,  J.  Opt.  Soc,  America.  44.  492, 

Glennon,  B.  M.,  and  Wiese,  W.  L.  1962,  Bibliography  on  Atomic  Transition 
Probabilities  (NBS  Mono.  No.  50  [Washington;  Government  Printing 
Office]) . 

Greenstein,  J.  L.  1948,  Ap .  J . .  107.  151. 

Hill,  A.  J.,  and  King,  R.  B.  1951,  J.  Opt.  Soc,  America.  41 .  315. 

Houtgast,  J.  1942,  The  Variations  in  the  Profiles  of  Strong  Franhofer 

Lines  along  a  Radius  of  the  Solar  Disc  (Utrecht:  Schotanus  6e  Sons). 

Hunger,  K.  1956,  Zs .  f .  Ap . .  39.  36. 

Johnson,  H.  L,,  and  Morgan,  W,  W.  1953,  Ap.  J..  117.  313. 

Keenan,  P.  C.,  and  Morgan,  W.  W.  1951,  Astrophysics,  ed.  J.  A.  Hynek 
(New  York:  McGraw-Hill  Book  Co.),  p.  12. 

King,  R.  B.  1942,  Ap.  J. .  9^,  78. 

King,  R.  B.,  and  King,  A.  S.  1938,  Ap .  J . .  87.  24. 

Menzel,  D.  H.  1936,  Ap .  J . .  84.  462. 

Minnaert,  M.,  and  Slob,  C.  1931,  Proc.  Amsterdam  Acad..  34.  Part  I,  542. 

Moore,  C.  E.  1945,  Contr.  Princeton  U.  Obs . .  No.  20, 


176 


177 


Schrpeder,  L.  W.  1958,  Ph,D.  Thesis,  Indiana  University,  unpublished. 
Str'djngren,  B.  1937,  Ap .  J . .  86.  1. 

Struve,  0.,  and  Elvey,  C.  T.  1934,  Ap .  J . .  79.  409. 

Swensson,  J.  1946,  Ap .  J . .  103.  207. 

Thackeray,  A.  1936,  Ap .  J . .  84.  433. 

Unsold,  A.  1955,  Physik  der  Sternatmospharen  (2d  ed.;  Berlin:  Julius 
Springer) . 

Wright,  K.  0.  1948,  Pub.  Dotn.  Ap.  Obs  ■  Victoria.  VIII,  No.  1. 

1962,  Astronomical  Techniques,  ed.  W.  A.  Hiltner  (Chicago:  The 
University  of  Chicago  Press),  p.  83. 

Wrubel,  M,  H.  1949,  Ap .  J..  109.  66. 

_ _ .  1950,  ibid..  111.  157. 

.  1954a,  ibid. .  119.  51. 

.  1954b,  Proceedings  of  the  Indiana  Conference  on  Stellar  Atmos¬ 

pheres  .  ed .  M.  H.  Wrubel  (1955),  p.  126. 

.  1956,  published  Aller  (1960)  above,  p.  199. 


VITA 


Hugh  Oscar  Peebles 
Candidate  for  the  Degree  of 
Doctor  of  Philosophy 


Thesis:  A  STUDY  OF  hXNE  INTENSITIES  IN  THE  SPECTRUM  OF  9  URSAE  MAJORIS 
Major  Field:  Physics 
Biographical : 

Personal  Data:  Born  in  Kountze,  Texas,  May  4,  1933,  the  son  of 
Hugh  0.  and  Gladys  L,  Peebles. 

Education:  Attended  grade  school  in  Kountze  and  San  Antonio, 

Taxas;  graduated  from  Kountze  High  School  in  1950;  received 
the  Bachelor  of  Science  degree  from  the  University  of  Texas, 
with  a  major  in  Physical  Education,  in  January,  1955;  attended 
Lamar  State  College  of  Technology,  Beaumont,  Texas,  for  two 
years,  beginning  in  1957;  attended  the  University  of  Texas 
during  the  summer,  1959;  received  the  Master  of  Science  degree 
from  the  Oklahoma  State  University,  with  a  major  in  Natural 
Science,  in  August,  1960;  completed  requirements  for  the 
Doctor  of  Philosophy  degree  in  May,  1964. 

Professional  experience:  Taught  in  the  Kountze  High  School, 
Kountze,  Texas,  from  September,  1955,  through  May,  195^. 


